


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations l. Thesis and Dissertation Collection, all items 


1974-12 


Stiffness of stepped bars 


Doyle, Dennis Michael 


Monterey, California. Naval Postgraduate School 
http://hdl.handle.net/10945/17253 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
) К D U DLE Y research materials and institutional publications created by the NPS community. 
FW : Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


Im KNOX appointed — and published — scholarly author. 

| | LIBRARY Dudley Knox Library / Naval Postgraduate School 

411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http: //wwwenps.edu/library 








NAVAL POSTGRADUATE SCHOOL 


Monterey, California 





STIFFNESS OF STEPPED BARS 
by 


верес тере 


mms 
--- -- αὶ ¿ —— w 


1 Doyle 


December 1974 


Thesis Advisor: R. E. Newton] 


Approved for public release; distribution unlimited. 





UNCLASSIFIED 


SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 


REPORT DOCUMENTATION PAGE 


. REPORT NUMBER 2. GOVT ACCESSION NO. 


4. TITLE (and Subtitle) 


Stiffness of Stepped Bars 








READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


3. RECIPIENT'S CATALOG NUMBER 







5. TYPE OF REPORT & PERIOD COVERED 
Master's Thesis; 
December 1974 


6. PERFORMING ORG. REPORT NUMBER 






. ΑΙ ΤΗΟΒ(ε) 8. CONTRACT OR GRANT NUMBER(e) 


Dennis Michael Doyle 


‚. PERFORMING ORGANIZATION NAME AND ADDRESS . ELEMENT, PROJECT, TASK 


ORK UNIT NUMBERS 
Naval Postgraduate School 
Monterey, California 93940 


11. CONTROLLING OFFICE NAME AND ADDRESS 
Naval Postgraduate School Decemberk ui 
Monterey, California 93940 AR 


4. MONITORING AGENCY NAME & ADORESS(II different {гот Controlling Office) | 15. SECURITY CLASS. (of thie report) 





Naval Postgraduate School lUnclassified 


EE [- a 


DISTRIBUTION STATEMENT (of thie Report) 


Approved for public release; distribution unlimited. 


DISTRIBUTION STATEMENT (of the abetract entered in Block 20, tf different from Report) 





16. SUPPLEMENTARY NOTES 


19. KEY WORDS (Continue on reverae eide if naceesery and identify by block number) 


otepped Bars 
Overall Stifiness 
Deflectjior c! Stepped Bars 


20. ABSTRACT (Continue on reveres oide if necessary and identify by block number) 


The stiffness of stepped bars is determined by calcu- 
lating deflections due to applied loads using the finite 
element method. Nine combinations of step height and 
fillet radius are studied ior each or the following cases: 








1. rectangular cross section under axial load; 
2. rectangular cross section under bending load; 
νιν. .::-: cress sGcticn uwaeer axial load; 
Е FORM | 
DD ¡jan za 1473 EDITION OF 1 NOV 65 15 OBSOLETE UNCLASSIFIED 


3 S/N 0102-014- 6601 —————————MM t 
(Page 1) а SECURITY CLASSIFICATION OF THIS PAGE (When Date Entered) 








UNCLASSIFIED Í 


€t CUHITY CLASSIFICATION OF THIS PAGE/ (When Όρια Entored: 


Block #20 continued 


4. circular cross section under torsional load. 
A correction parameter, actually a fictitious axial dis- 
placement of the step, is derived. It makes possible 
accurate calculation of deflections of stepped bars 
using elementary formulas. 


"M 173 ACK 
BEES “5 Es ; UNCLASSIFIED 
S/N 0102-014-6601 SECURITY CLASSIFICATION OF THIS PAGE(When Data Entered) 





Stiffness of Stepped Bars 
by 


Dennis Michael ‚Doyle 
Lieutenant Commander, United States Navy 
B. S., United States Naval Academy, 1966 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN MECHANICAL ENGINEERING 


from the 


NAVAL POSTGRADUATE SCHOOL 
December 1974 














DUDLEY KNOX LIBRARY 
NAVAL POSTGAEUATE SCHOOL 
MONTEREY, CALIFORNIA 93940 


ABSTRACT 


The stiffness of stepped bars is determined by calcu- 
lating deflections due to applied loads using the finite 
element method. Nine combinations of step height and fillet 
radius are studied for each of the following cases: 


1. rectangular cross section under axial load; 


2. rectangular cross section under bending load; 
3. circular cross section under axial load; 
4. circular cross section under torsional load. 


A correction parameter, actually a fictitious axial dis- 
placement of the step, is derived. It makes possible ac- 
curate calculation of deflections of stepped bars using 


elementary formulas. 
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There are many mechanical and Ет applications 
of bars and beams which have abrupt changes in cross-sec- 
tional area. In the design of these members their deflec- 
tions under expected loads must be considered. If accuracy 
is not essential, the deflections can be determined by 
relatively basic calculations. These calculations involve 
formulas based on the simple stress distributions found in 
constant section members. They do not take into considera- 
tion the reduced stiffness due to stress flow irregularities 
in the vicinity of the step. A stress concentration is said 
to exist at this location. 

Through the numerical analysis technique called the 
finite element method, the deflection of stepped bars or 
beams can be accurately determined. This method is sensi- 
tive to the effect of the stress concentration. 

The purpose of this study was to analyze the overall 
Sctiliness of stepped bars. An additional goal was to de- 
termine a correction parameter which could be used in de- 
sign to compensate for the reduced stiffness due to the 
stress concentration. Finite element analyses were performed 
for stepped bars of various dimensions under axial, bending, 
and torsional loads.: The resulting deflections were deter- 
mined, converted to a common correction parameter, and 


compared. 


ł0 





The recent emergence of the finite element method along 
with the digital computer has made feasible the previously 
prohibitive calculations necessary for stiffness analysis. 
The only prior study of stiffness of stepped bars that this 
author could identify was accomplished by F. P. Porter Sa 
His study of torsional vibration of irregular shafts in- 
cluded stiffness calculations based on the mathematical 
theory of elasticity using an approximate graphical 


"construction. 


sumbers in brackets refer to the List of References, 
ES O... 
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πι. -ΙΠΡΕΜΕΝΤΑΚΥ STRESS RELATIONSHIPS FOR THE THREE LOADING 

EN ud c °° °°» 
A. BASIC ASSUMPTIONS 

For each of the three loading cases (axial, bending, and 
torsional), the stepped bar was assumed to be perfectly 
elastic, homogeneous, isotropic, and to exhibit a linear 
stress-strain relationship. In addition, it was assumed 
that the bar was initially straight and all loads were ap- 


plied smoothly to avoid the effect of impact loading. 


B. THE AXIAL LOAD. CASE 
A bar subjected to equal and opposite tensile loads will 


deform according to the equation 


= е en 
where ó is the elongation in length £, P is the axial ten- 
sion, A is the cross-sectional area, and E is the modulus 
of elasticity of the bar material. This equation is valid 
. provided that the loading is appiied uniformly across both 
ends of the bar and that the bar is of constant cross-section. 
For a bar of non-uniform cross section, the elongation 

is estimated by dividing the bar into segments of constant 
cross section and adding the segment elongations. This pro- 
cedure does not take into account the reduced stiffness due 


. to the stress concentration located at the change of cross 


Lon. 
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Figure 1 is a longitudinal section of the bar studied. 


The estimated elongation of the bar is therefore 


_ Р 
6 =-F[ Z + 1. - (2) 


А, апа A, are the respective cross-sectional areas. 
The axial load case was analyzed for both the plane 
stress state and the axisymmetric stress state. In the 


plane stress analysis, the bar was considered to have unit 


ia 1. —а 


г 


Figure 1. Geometry of stepped bar. 


thickness, therefore the cross-sectional area equaled the 
height of the bar. In the axisymmetric stress analysis the 


pornas considered to be of circular cross section. 


Fair SANAA AA Oe 
DET BENDING CASE 


o 


For bars subjected to flexure, or pure bending, a mea- 
sure of the beam flexibility is the angle through which the 
ends of the bar rotate. The total angle can be determined 


from the equation 


mE 

Hilo 
— 
CO 
м2 


ф = 


1. 





where $ is the sum of angles Ó; and Ф, as described in Fig. 
2. M is the bending moment and I is the moment of inertia 


BI the horizontal cCentroidal axis of the Gcross—-sectional 


Figure 2. Bending angles. 


area. Equation (3) is derived from the following equation 
which relates the bending moment around the neutral axis of 


an elastic beam to the curvature, kK, of the elastic curve. 


M 
Bee 


K = 
This theory is valid provided: 
1. Plane sections remain plane during bending. 
2. The bar is ori constant cross section. 
3. Bending is such that shear strains are negligible. 
As in the axial load case, the stepped bar can be con- 


Sidered in two sections, each of constant cross section. 


The estimate for the total bending angle is therefore 


> 
πι 
> 


2 


2 


+ 


ПЕ (4) 


X» 

1 
mz 
m 
| 





H 
= 
H 
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The bar analyzed was of rectangular cross section and 
considered to have unit thickness. The bending moment was 
applied by a linearly varying normal stress across both 
ends of the bar exerting tension on one side of the neutral 
surface and compression on the other side as shown in Figure 


з. 





Figure 3.  Bending stress distribution. 


Bee LHe TORSION CASE 
The equation for angular deflection of a uniform, round, 


crossesection bar subjected to torsion is 
ο. 
0 = s `: (5) 


ne angular deflection, T 3s the torque, J is the polar 
moment of inertia of the cross-sectional area, and G is the 
modulus of elasticity in shear. It is assumed that plane 
sections perpendicular to the axis of the bar remain plane 
and that a radial line remains straight when the bar is 
twisted. 

In this study, one end of the bar was held fixed anda 


zoraue was applied to the other end such that the shear 
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stress varied linearly from zero.at the centerline to a 


maximum at the outer radius as shown in Fig. 4. 


Figure 4. Torsional stress applied to circular cross section. 


As before, the bar was considered in two sections of con- 
stant cross section and the total angle of twist was the sum 
of the angles for each section. Therefore the estimate for 


the entire bar is 


+— 1]. (6) 


16 





I. ο Μο ο ο Ες σ КТБОНМ 


It is well documented" rm Section of a bar of 


non-uniform cross section is slightly greater than would be 


indicated by basic equations (2), (4), and (6). 


his dez 


creased stiffness is caused by the stress intensification 


at the point of discontinuity. 


One way of accounting for this decreased stiffness is 


teomconsi der the effect of the stress intensification on the 


bar as a whole to be the same as the effect of moving the 


location of the discontinuity. Increasing the length of 


the smaller portion of the bar and decreasing the length of 


the larger portion by the same amount would achieve the same 


Section ts stitfnesso 


This apparent change in the location of the 


bar discon- 


tinuity was calculated and used as.a measure to compare 
stiffness of bars of different dimensions. The size of 
this change was defined as the variable "A" and is illus- 
reden Fig. 5. 
ο 
222 
k— - — δη ----------- {ο -——— e 
Figure 5. Geometric interpretation of stiffness correction 


parameter. 


do 





To calculate A, equations (2), (4), and (6) were altered as 














follows: 
JA k +A - 
P ab 2 
ó--l + ] (7) 
E Ay A, 
~ А £ +A 
φ-ξῖ--- +<] (8) 
n 2 
2, -A КРА 
en +) (9) 
n 2 


In each loading case, the deflections required to cal- 
culate 6, $, and 0 were determined by digital computer finite 
element stress analysis. Knowing the load, modulus of elas- 
ticity, and physical dimensions of the bar, the parameter A 
was determined. The A parameter is common to equations (7), 
(8), and (9) and therefore is a means of comparing stiffness 
corrections for the three types of loading. In addition, it 
15 independent of the load applied because the magnitude of 
the deflection and the load applied are directly proportional. 

Since the A parameter is a means of correcting for the 


deflection of a stenmned bar under load, it was defined as 
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IV. SOFTWARE 


A. BACKGROUND 

Before the stiffness correction parameter could be de- 
termined, the deflections of the bar under each of the three 
loading conditions had to be determined. Although photo- 
elastic analysis and direct strain measurement are suitable 
for determining the effect of stress concentration, the ideal 
method for overall stiffness evaluation is the finite ele- 
ment stress analysis. 

In very brief and general terms, the finite element 
method is an advanced numerical method that has been made 
possible by the high speed digital computer. When applied 
to structures it can be used to perform two or three di- 
mensional stress analyses. The structure under study is 
divided into subdivisions called finite elements. Adjacent 
elements are connected along the entire inter-element 
boundary. Specific points at regular intervals along the 
element boundaries are designatea and derined as nodal 
points or nodes. When a load is applied to the structure, 
the displacements of each of these nodal points are calcu- 
lated and used to determine the stresses. 

For each of the problems, maximum use of existing com- 
puter programs was made. For the plane stress analysis of 
the axial and bending load cases, a program titled PLISOP 
[2] was used. It was developed at the Naval Postgraduate 


School, Monterey, California and has been in use for several 
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years. For the torsional loading case, a program titled 
TORT 2 [3] was used. It was developed. at the University of 
Wales, Swansea, and had been disiaiseu locally but had not 
been used or completely debugged. In order to study the 
fourth problem, the axially loaded bar of circular cross 
section, the capability of PLISOP had to be expanded to 


perform an axisymmetric stress analysis. 


B. ΜΕΡΗ GENERATION 

The network of lines describing the division of a struc- 
ture into finite elements is referred to as a mesh. Nor- 
mally a finer mesh more accurately describes the physical 
structure under study and results in a more accurate stress 
analysis. Unfortunately, increasing the number of subdi- 
visions or elements greatly increases the computational ef- 
fort involved in the problem solution. As a result, the 
maximum capabilities of most er purpose digital com- 
puters are rapidly reached as meshes are refined. 

In many stress analyses, the area of interest is only a 
small portion of the structure. In these cases the mesh 
can be very fine in the area of interest and much coarser 
for the rest of the structure. This minimizes the total 
number of elements yet achieves accuracy where it is desired. 
In this study, however, the analysis of the structure as a 
whole was of interest. Therefore it was desirable to obtain 
a mesh as uniformly refined as possible within computer pro- 
gram limitations. The only practical means of maximizing 


accuracy yet minimizing the total number of elements was to 
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take advantage of symmetry. This was done by dividing the 
bar in half along the longitudinal axis. 

It was found satisfactory to take the lengths E he ος 
the two bar segments equal to the larger height dimension 
D. This provided a nearly undisturbed stress distribution 
at both ends of the bar. Thus the stress intensification 
resulting from the abrupt change of cross section is negli- 
gible at the ends of the bar. 

Based on these considerations, the stepped bar was 
represented as shown in Fig. 6. In each loading case the 
ratio of d/D and the fillet radius were varied. These vari- 
ations necessitated minor changes in the mesh representation, 


however, this general discussion will refer to the arrange- 


ment of Fig. 6. 





Figure 6. Finite element mesh. 


For this basic geometry, the mesh used consisted of 
twenty-four isoparametric quadrilateral elements. The con- 
cept of the isoparametric element requires lengthy ezplana- 
tion and is not considered essential to this development. 
The reader is referred to chapter eight of Ref. 4, for the 


background of the isoparametric element. Quadrilateral 
Zu 





elements were used as opposed ος» elements. 

Both of the computer programs used isoparametric quadri- 
lateral elements. Such elements can be linear, quadratic, 
or cubic. In this study, the only curved portion of the 
mesh, the fillet, was adequately described by quadratic 
elements. 

The quadratic, isoparametric, quadrilateral element has 
four corner nodal points, or nodes, and four mid-side nodes. 
The coordinates of each of these nodes are part of the in- 
put data for both the PLISOP and the TORT 2 programs. The 
task of providing this coordinate data as well as data to 
describe how the elements fit together is a considerable 
one. Fortunately, another computer program was developed 
to perform just this task. This program, titled PLIMEG [5], 
generates a mesh as specified. It requires only enough geo- 
metric data to describe the overall area to be discretized. 
The directions for use of PLIMEG are contained in the program. 

The twenty-four element mesh was chosen as an adequate 
representation of the region. A later convergence study, 
which used a finer mesh, confirmed the adequacy of the 
twenty-four element mesh for this primarily comparative 
study. 

Illustrations of the meshes used for the various geo- 


metries are shovn in Appendix.A. 


ES COMPUTER PROGRAM EMPLOYED 
l. The PLISOP program performs a plane stress or a 


plane strain analysis by the finite element method. Detailed 
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directions for its use are contained in the program. The 
bulk of the input data" cetmebeosperopone? Бу the mesh generat- 
ing program, PLIMEG, previously mentivired. The axial and 
bending effects were achieved by the loads applied and by 
the boundary conditions imposed. 

The load vectors must be derived externally and are ap- 
plied by specifically defining the load and its direction 
at particular nodal points in the mesh.  Derivation of the 
force vectors for both axial and bending loading is shown 
in Appendix B, sections 1 and 2, respectively. 

The boundary conditions necessary for the axial loading 
of the bar half-section shown in Fig. 6 consisted of con- 
straining. each node along the bottom surface (bar neutral 
surface) to freedom of movement in the longitudinal direc- 
Ν᾿. a addition. one of the corner nodes along this 
surface was pinned to prevent longitudinal motion as well. 
For the bending load, the nodes along the bottom surface 
were free to move in the vertical direction only. The two 
corner nodes along this surface were pinned to allow only 
rotational displacement at these two points. 

2. The TORT 2 program performs a stress analysis of 
circular sectioned solids with varying diameter using the 
finite element method. The directions for its use are 
contained in a separate booklet accompanying the program. 
The same meshes were used to represent the stepped bar 
except the origin of the coordinate system had to be at 
ου σος right corner of the bar. The origin for the mesh 


used previously for the PLISOP program was at the lower 
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left corner. A relatively πρ τας routine can be 
used to renumber the nodal point data provided by the PLIMEG 
program, but the connectivity data was revised by hand. 

Loading is applied by specifically assigning values of 
shear stress at particular nodal points. 

3. To study a bar of circular cross section under axial 
load, an axisymmetric stress analysis capability was re- 
quired. The most direct way of obtaining this capability 
was to alter the PLISOP program. Due to the symmetry of a 
body of revolution, the problem remains mathematically two 
dimensional. The principal difference between the plane 
stress and the axisymmetric stress problems is the intro- 
duction of a fourth component of stress, the "hoop" stress. 
The details of the development of the axisymmetric capability 
are included in Appendix C. 

Since the PLISOP program is essentially unchanged by 
this added capability, the directions for its use are the 
same. Appendix D contains notes on the axisymmetric capa- 
bility intended to amplify the directions contained in the 
program. Appendix E contains a listing of the PLISOP pro- 
gram in its revised form. 


tress analyses of the axial and bending 


- 


A ә ү ET 
ае Salome 


0 


load cases, the load vectors must also be derived externally 
for the axisymmetric stress analysis. This derivation is 


included in Appendix B, section 3. 
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V. RESULTS 


eee ERE SENTATION OF RESULTS 

Bon each of the loading cases, nine variations of the 
bar dimensions were analyzed. The dimensions and results 
were nondimensionalized by dividing each by the dimension, 
D, which was held constant throughout the study. The 
values of the stiffness correction parameter, A/D, for each 
of the four problems are presented in tabular form, Tables 
ἵνα graphically, Files. 7 - 10. 

From the graphical presentations it is evident that the 
largest correction parameters are required when there is no 
fillet. As the fillet radius is increased, A/D decreases. 
In addition, as the dimension, d, is varied with r held 
constant, the largest correction is required for the ratio 
of d/D equal to one-half. 

The negative values of A/D indicate the bar is stiffer 
than equations (2), (4), and (6) would show. This is ex- 
plained by the ract that these equations, as applied, do 
not take into consideration the increased cross-sectional 
area resulting from the fillet. The fillet area was not 
included in the equations for calculating the correction in 
order to keep the calculations simple and uniform for all 


ases. 


EM CDUPARISON OF THE FOUR CASES 
In order to compare the four cases, a plot of A/D versus 


r/D was prepared for each of the three ratios of d/D. Figures 
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Table I 


A/D for Axial Load 
(rectangular cross section) 





Table III 


A/D for Bending 
(rectangular cross section) 


СТОЯ | | | 
A 


r/D "3. PEN " 
ΠΗ; Ι΄ 


Dr 10.12 [0.11 | 
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Table II 


ADO Aia Ll Load 
(circular cross section) 





Table IV 


AD foc Torsion 
(circular cross section) 


Na | | | | 
e N28 [0.50 [0.76 | 
| | 


ο 0 020371006 
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-0.20 


Figure 7. Stiffness correction for axial load (rectangular 
σσ σε τοι). 
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Figure 9. Stiffness correction for bending (rectangular 
cross section). 


29 





ci" 
N 
O 
e 
O 


0.15 


| = 0.125 
0,10 


0.05 2 DC 2 Г 5 = 0.250 


À O 
2 0.0 


η 


-0.05 


-0.10 Е 


-0,15 


-0. 20 


- Figure 10. Stiffness correction for torsion (circular 
cross section). 
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mes 12, and 13 show the effect m stiffness of the three 
different types of loading and two types of cross section. 
The most severe type of loading is clearly axial as indi- 
cated by the larger values of A/D. I. M. Allison [6] also 
found the axial loading to be the most severe of the three. 
His conclusions, however, were based on stress concentra- 
on actors. 

Comparison of the two types of cross section subjected 
to axial load indicates the circular cross section is 
stiffer than the rectangular cross section for a particular 


po of d/D. 


C. CONVERGENCE AND UNCERTAINTY 

The normal methods for proving convergence of a finite 
element solution are by further refining the mesh or by 
using higher order elements. Both of these methods were 
investigated. For the basic bar case (D=4, d=2, r=0) under 
axial load, the convergence test results are presented in 
Pig- 1}. 

These results show the improvement in accuracy gained 
by uSing the quadratic rather than the linear elements. Al- 
though it can be assumed that cubic elements would give 
further improvement, other considerations prevailed. Since 
it was desired to cover a wide range of problems, the addi- 
tional computational effort required for cubic elements was 
deemed an unnecessary luxury. Also, since the geometry of 
the bar was relatively simple, consisting of plane or cylin- 


dric surfaces, cubic elements would not provide the same 
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Figure 14. A/D vs.l/(no. of nodes). 


percentage improvement that the quadratic elements did. 
In order to check the accuracy achieved by mesh refine- 
ment with quadratic elements, the PLISOP program limits were 


altered to accommodate a ninety-six quadratic element mesh. 


Γ t] 
0ч 


σος 19 in aces this mesh Las come very close to con- 
vergence. Based on these results, the uncertainty in the 
values of A/D in Tables I - IV is considered to be the dif- 
ference between the ninety-six element result and the twenty- 
four element result, rounded to the second decimal place. 

The uncertainty is therefore +.02. 

An additional test problem was run to verify that the 


geometry studied, with dimensions 2, and 2 equal to the 


2 
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larger height dimension, was adequate to ensure that the 
stress intensification at the discontinuity did not signifi- 
cantly affect the stress distribution at the ends of the bar. 
To check this aspect, the length dimensions E. and Le, were. 
increased to 1.5 times the larger height dimension. The 
resulting A/D parameter was within .001 of the corresponding 
case with the shorter length, thus confirming the sufficiency 
of the ratios 2, /D=2,/D = 1.0. 

The one geometry and loading case which can be compared 
with Porter's work agrees very favorably. For torsion of a 
Ber with d/D = 0.75 and r/D = 0.0, thevstitiness Correction 


parameters agree to the third decimal place. 
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VI. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

The stiffness correction parameters obtained are useful 
and applicable in the design of stepped bars. The tables 
and curves included in Section V can be used with confidence 
to the first decimal place of A/D within the range of physi- 
cal bar dimensions studied. 

The results also yield important comparative data for 
stepped bars of various dimensions under the three types of 
loading. The following comparative observations were made: 

1. The value of A/D for stepped bars subjected to axial 
load is larger than the corresponding value for bending or 
torsion. 

2. The stepped bar with ratio’ d/D equal to one-half 
yields a larger value of A/D than the other ratlos studied. 

3. Fillets reduce stress intensification at the loca- 
tion of the step and therefore increase the stiffness of 
the bar. 

In addition, it was observed that the localized stress 
intensification had a negligible effect on the cross-sectional 
stress distribution at a distance equal to the height of the 
Bar. This illustrates the St. Venant principle of rapid 


dissipation of localized stresses. 


B. RECOMMENDATIONS 
It would be of interest to isolate one loading case and 


through use of more refined meshes and cubic elements achieve 
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correction parameters accurate to the third decimal place. 
In addition, a great deal more.caomnav2tive information 
could be gathered with a three dumeisiconal finite element 
solution program. Both of these would require altering 
existing computer programs and would greatly increase the 
computer core storage space required for each run in this 


study. 
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APPENDIX A: 


FINITE ELEMENT MESHES EMPLOYED 
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APPENDIX B: DERIVATION OF EQUIVALENT NODAL FORCES 


1. Plane Stress Analysis of the Axially Loaded Bar 

Boundary stresses and distributed loads are converted 
to equivalent nodal forces by the principle of virtual work. 
For the case of external stress o applied along the boundary 


of an element, the nodal forces are developed as follows: 


ôW = [δυ.1΄Γε.15 = f o*óu. d(area) 


substituting 


TEN SOL JI ENSIS 
1 1 
[£,J% = $ o[N,Td(area) (10) 
where 
6W = increment of virtual work 
f. = nodal forces 


o = applied stress 


Su, - incremental nodal displacements 
δι = the virtual displacement of a typical point 
on the element surface 
N; = nodal shape functions 
e = superscript, denotes element contribution. 


The quacratic, isoparametric, quadrilateral element 
pictured below was used throughout this study. The shape 
functions are expressed in terms of the local, normalized 
coordinates (&,n). For loading as shown, the nodal force 


equations become: 
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Figure 24. Quadratic, isoparametric, quadrilateral element. 
ae БЯ N, 
= E 
-] 
f N. 
where 
h = height of the element 


d(armea) = dy = (h/2)dn for bar of unit thickness. 
Substituting the shape functions and integrating, the result 


for c equal constant is: 


m А 
2 Eje 

e E oh 

F 

ge y 


Where nodes are shared by two adjacent elements, the 
nodal forces are summed. 
2. Plane Stress Analysis of a Bar Under Bending Load 
The equivalent nodal forces for the linearly varying 


Stress were derived in the same manner as the axially loaded 
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case except for the applied stress. The linearly varying 


stress shown on the element surface, shown below, can be 


73 
7:0 To) 
a= σι 
Figure 25. Elemental linearly varying stress distribution. 
represented by the equation 
с = с, =" τι 1 ). UIS 


T РЕ FN. 
- l-n 1+n y,h, 
f, = / С 5 PEt σαί > ) > N, an 
= | 
fj η 
τη xe 
f - 2 fog * 2g (13) 
2 6 1 σα | 
fi τ 


Applying this formula to one end of a particular bar 


Spanned by two elements as shown in Fig. 26 and letting 0470, 


ey 25g 2 
then ση - σα 2399/2 and 01 


Therefore, the nodal forces along this surface, num- 


= 0. 


bered from top to bottom are: 
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superscript = element 





number 
subscript = element 

node number 

Figure 26. Global bending stress distribution. 

1 
Р, f. 1 
Ε ү 3 
2 1 2 ho 

Pa = т + f3 "cS l (14) 

Fy x I 

He Les ] o | 


3. Axisymmetric Stress Analysis of the Axially Loaded Bar 
In this case, the region associated with an element is 
the sector of an annulus with a one radian central angle. 
Figure 27 illustrates the coordinates used with the axisym- 
metric element and the node numbering used in development 


below. 





Figure 27. Axısymmetric element. 


ου 





Beginning 


for a uniform 


with equation (10), the equivalent nodal forces 


external stress are developed as follows: 





[Ε.1- Ξ [o-[N, Jd(area) (10) 
d(area) = r (l radian) dr = г dr 
r= суйу + rN, T rN, 
= t t t 
dr = (түт + rN, + rN, )dn 
where N! = 3N./9n. Then, 
1 1 
r. =- r 
dr = poe апе: 
2 
Therefore, 
Er cr 
d(area) - > + rN, + r. N.)dn 
Substituting back into equation (10), 
2 
г, f (NJ) dn FN¿N¿dn fN N.dn rj 
r. , 
f, SC pe S(N,) dn SN Naan ry 
2 
f. sym. f CN4) dn r, 
integrating the shape functions si DET leaves 
І y 2 -l ES 
rr i 
f, acus 30 215 2 r, (159 
f, -1 E y r, 
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ΙΙ DEVELOPMENT OF THE AXISYMMETRIC STRESS CAPABILITY 


For an axisymmetrically loaded bar, displacements are 
limited to the axial and radial directions. In this study 
these displacements were identified with the letters u and 
v, respectively. Cylindrical coordinates are applicable 
and are defined in Fig. 28. It can be seen that radial dis- 


placements introduce a circumferential component of strain 


r 





Figure 28 Cylindrical coordinates. 


as well. as a radial component. For problems of this type, 


the general strain-displacement relationships, as developed 





in chapter five of Ref. 4, are 
r € ο] Л ο 1 
Z 
E | dv/or 
- : (16) 
a v/r 





ου οσα. Αμ даг 


Bes s К = 
bel | 


The elasticity matrix which links the strains to the stresses 


1s derived from the following stress-strain relationships: 


Lo, - v (o, + σρ) 


N 
ΓΙ p 


io. = а саи 
E Z 8 


о З 





С WI 
Eg = =F [σρ - v (o, + аз) 
Yoz x Τα’ 6 


The resulting elasticity matrix [DJ] is 


1 v/(1-v) v/(l-v) 0 


j v/ C1-v) 0 
N E(1-v) 
[D] = (1*v) (1-2v) sym. 1 0 m 


(1-2v)/2(1-v) 


The finite element formulation of the element stiffness 
matrix is carried out as in the plane stress problem, ex- 
cept the integration over the volume of the element must 
take into consideration the fact that the element thickness 
varies with radius. The element stiffness matrix [k] can 


therefore be expressed by the equation, 


кай К а лу л 2 (18) 


where [B] is the matrix derived from the strain-displacement 


relationships and defined such that 


e 


[e] 5 ἱδ]ία] 


where [o]^ is the element displacement vector. 

The primary changes to the PLISOP program were to in- 
corporate the new [B], [D], and [k] matrices, as defined 
above. The element shape functions also had to be added, 
since only their derivatives were included PLISOP. 

For further detail on axisymmetric stress analysis by 
the finite element method, the reader is referred to chapter 


five of Ref. Ч. 
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EXUSENDIX D: NOTES ON AXISYMMETRIC STRESS ANALYSIS CAPABILITY 
OF PLISOP 

l. z and r denote respectively the axial and radial 
coordinates and u and v the corresponding displacements. 
ο του z Z x and р = у.) 

2. The volume of material associated with an element 
has a thickness equal to the arc subtended by a one radian 
angle vice the body of revolution usually discussed in most 
texts. It turns out the only difference is a factor of 2m. 

3. Since the axisymmetric analysis involves division 
by r, there cannot be any nodal points with coordinate r=0. 
Fopr nodal points along the centerline, substituting p-1x10^" 
suffices to prevent the solution from "blowing up," although 
the accuracy of stress values near these nodal points is ad- 


versely affected. 
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APPENDIX E: COMPUTER PROGRAM 


EOOOOO!UnooooooooooooocowuwoooOkt«QOooceoOoooooooooooooooo 
CUA ЧОМ МО СО — CN σας ων οσο ο συ ον ο CDS I fS- 000 Оч СУ XT UO O OO — CN 
OO OO OO O O O c c rA c mA mH mA == С\] С\З с\у СМ AICI AIA СМ СМ 0e 00 00 e СС СС) «T «Xr E Jr Xr sts ELA LALA 
C O OO O O O O C C C O C C C C O CO COCO O CO O O OCO CO O CO O O O OC O O COCO O OO O O O OO S O 
G OOS) 5 OO CO OO CO OOQOƏCOCOƏOOOCOƏOQOƏOCOƏOGCOƏCOQOCOOCOQOCOƏCOGCOGOOGOCOOƏOOƏOƏOƏOƏOO 
22 CZ C LL ZALLA LALALA ZOO ZE E ZA 22 Z ZZ Z ZZ xx= A АСР A AZ 
<( < Z= — Z <1I X` <( <( <( <( <( <(<( <( <( <í <( <( <( <( — > <( <( <( <( <( <( <í <( <( <( <( <( <í <( <( <( <( <( <( <( <( <í <( <( <( 


DOAAIVDO OO OO OO OOVIOOOODIOOO VO OOO OO OD OVOOOVOOOOOOOOO 
+ 


х e 
X X x c < 
X om jt 6) O цы С 
+ +++ =t [== LJ ο 
T* > M. 3 Dom = O σας 
* | — ο =) ш O ш us cor zx 
+ = I ir Dem Z =- YN > «Ш> ш 
+ >< e > I ue 7 = „ou 3 
X wo 9 +£ = = ae =a INFOZO ш O © m 
JE »-ν- Ox ш ш OX OQ шши <A m azz o 
FRA UJ 3t ч JJ мис OXUXZO с м< of 
+ IWW O 3 D da e Met Woetlwwii O ολ = a 
+ = + O O e ZO ruLror o bt LLY 
X ZD LLJ3£ с X ο Оц ZIO mk A -Aqa Ες 
y < Z -3t ο. O. FE e — eR UU Z = UJ 
+ -- »-Ἑ 22 UN OOZ ew cx wan a 
Hee eD OF oc ο wo — 2 U3UJ oor ш E «Ie 
+000 A+ HJ) us CAO LUT ο oO. 2? и) 
+ «OT 546 G. O. L Z CeO Ze > < e 
+ ZILIO > «al E O "kh ZUTA υ) vd» ax 
ei 2) κ. Cc) CÓ) шо -— Quz uwucdrz Aa = <q 
κ σιυιο Ον c A u O = atk oO c m= WO0 OO 
К ССИ X <“ < N ο) οι) ZWO Q < ZZ o zZ 
4e FT w+ O Ω Q em EI CY LU j O pu Y «Y 
+оо OX =ч PS ώσασωαωώωιωα. Ζω CO ZU = 
+ Z Ox ш ш πια κ) ο NWT rm “STO < 
4 uJu0c „+ Z Z worn JZU m e CX zx = 
к 25-0] + O Mee: . Sansa ш = шшс С > O 
%+ < 007) >- nxo0O Zw τι ος тос 3 =ч 
4. J (46 θ TI ON) OM OULD bp pon O) but Un 
KAUW + x δι NY [ml O J = >. 
+ e UJ OX «1% ewel Tie Quo — = < cN < 
% NWUXWALI + SOX Шии Ооо 5 N MTSGu- ο a. 
3£ C) FO <( t— C š чї eo БШО e zZ X CC ο x 
um DOG —IZ= JZ ZG —OWNEY =ч COL e e Um 
3 K Z Ç) — < + UND mm <A U) OJO «ru e СО. x иц = 
3r Er + = =QO)NuUOA dc. O o acc + „U 
FI ZI>* = 21002 2и «Ο «2 9 --. Όι = < > 
+ LW ke + = WOFWOZ< naLe NA uJuJ O O = Ww 
+ ШШ oce H < UY XUZ E (QG*-OZX Eu e IIlIuu und = a 4 
εστω. O = ш mel OA TDL Сс) 0 ehe I. The > ë шш 
3 «E .IEst N O ! Ош ци πω" ш + OÖ = SI 
zz U Qs w ~ WM ο ασ OTO- (ul e αμ L zu ~ = LJ 
yt O.> <I + ~ ου. «ους πο ΟΖ ο OOUU Ох eu I 
x MW -»0 + LL >< UL O Ou. Z L.L = 4 ООО м U D ц. bk 
Aa mZ ay Ооо ш =O A € LL UJ mC) t^ JO 
$ NLO% = хохо co = ekke «o E uL 
HE YN» 4ς XA CWSU X C LL U. u. weet Z Z 5 ( < > C О 
+ = Annu mox W eL O D O U þe нан 2 ἡ OZ Zw O ZO Z 
+ C >> «1020 ә Бл о. Рн ЧИИ ЧОРО ЕС TT os oe geet) Ἢ E O 
3 O —I<T COLLI i£ suus D > ni ος ος Gar E VU) cus wu) 
FU IX mt O у и τε leet ot pet E =Zr cr O. Бш uy «(uuu 
3: cé Z < ¿— % ш τουσ ο -α uj QOO в MT ws мм») 
+u<aawons O ш ш < LU hm uU) We OOOO Ne KF. > C x 
a ОСО + Lu M J mm OW ry O>Z ZO AZO вә с ‚ш уц. о) 
4 D ш πο. [in le DJ. © ML C) < шоњ тна 
A + Z= AU) = OC (5b Eu < <t O ZO G S T 
el) DX% «νεο 222 ge μα. II Km ee X U. επ OO 
UL X OD sé Y eO or U^ x en c 
Home Wyse o < + O< o. < 
«Όλος: x e. > ος -- 2» 
кми чок qa <{ 26, 
++ OQ.. C) C C) >= KH i: = ο) < 
A OLUN A UJ ONONO Nono OWOWMOW м ΟΟΟΟΟ 
+ «ХЦ ч Н J ON ROUND OO PEL bt A UN 5 СС eet NITY 
Ἡ μαι БО» D Uu O| | ! Í || tree aof tl ry i uu) titel 
$ TUZAX AS a -2 OQ- ONOHO чо-чо rl O) 0) ч ee O —— ——— 
ἐξ — >` — Q. — + Z m= ος) 4 4 ЈСО MMS OD чеччи τς) nm. 

— — a u) A. 

E 


QOQ(QQOGOQOQOGOOCOOCCOCOCOC ODOODOOOOOODODOO DO DO DIOOOD OO DODOOODOOOOD 


96 





Q CO CO C O O O O C O O CO O CO O O O O O C O CO CO O CO CO CO O O CO CO OCO OOO CO CO OO CO O O OCO OO 
(o d try NO — OO СТО іф омс G" O — (NI r QF LO T.O I— OO DON cO +f try O — OO ON O = (NJ co +f О М со 
LAVAL LA LA UN NO WOO WO 0 0 0 Оо Pe Pe Da De Due De De De Pe ICH CO O cO O O O O Q QA G AQ A ON OA OA OA CA O, 
cid le C, a O O O C O C O CO O OO CO CO CO CO COCO CO OCO CO CO C O CO CO CO COCO OO COCO CO CO OCO OCO OOO O 
C OO CO CO CO O O O CO O O Co C O CO 090 95 995 9 99 995999599909 9 9 99959 9 959 9 OO 9 9959 
“>> 122 22 Z 2222221212 222 г> г УУ г = > οο 
<L <I <1 <I <T <I <I <1 «T«r «T «T«tC«[ «T «1 «C «r «T«r «I«ep[«I«[«1 «T «[Q«t1 «t1 «T «t «1«1 «1 «1 «1«1«1«1275 =— Z = X 
ENUDUODOOOQOQOOOOOQOQOOOQ(QOOOQOOOOO0OU0O0O pDDOOO0O0O QOoOoOo0oU0occococo 
OO 9 LL. 
cc bL ° © © 
шш Ç) hb мМ un 
NN ш U | Ou Luz 
=. ox ш ο < С < 
”' ш UU) CC X min <> N O 
LL y „4а Ц) СӘ кы e Lu DD № “ч 
= со J O ошо» > _ © (3 
co e oe > ZOU «k Lu PO =ч Lu 
© сс с Б Шш O»xiucoulzm = NO = =< 
ας ш а QoccuJh- (с сс C ceo <I μα) 
а. eo Wet Thr! >< (O G Z= O Son -- > 
= m AX wm Ole ZU] < xX о О oO 
ος κ Ü чы) Б ССО < С. шо О иш 
шш = a. OOZ = QA ти I Eo) 
а. ze er OMe wT α μα ы = 
Z Y >- 2<ц! YON Ce | U) = YN ши. 
2 О O zZz Ὁ ο ιο. a hum OT 
ο m e cc a ZZOa O OZ ~ Jm WY) Dæ ш = 
ес U PZr < mest OF A μπι στο oc Mm m o (Y 
<t Nk Zu Ὁ ze, z Z= e — eO б и < > O 
Oo M ZUG e *s DQO" OJO „Jr (©) > =) ei, ts LL) 
Z UJ Z UJ e= ano ο. ο uio! «t .J Cul C >< - 
YU. Ju 0 YAN WO (CC Z TOO vo шт «<= — — pao (0) 
aJ OW Jus с ---- e° T «КШ у ШШ pb tN ES Yr c ο Z= z 
шш Ό-ω = б OOO OO kea I eam Z 2 ᾽-- ыз (9 CO O UJ << 
= =U eZ Z NM δι ος 2ο) οι è ә <I = = = 
oO —ZUJ a Om > UJ ZAN a ZD =Z VY) Quo 
CL K — Z — G ANNOY te >j O O m = OOO eo ne O — O < 
uu ZO mt к-г ЛС XAO ч а. Z Ug 
Ore OO >< Te ee NTO C LLI LL | = 
SADA. .ἠ ite <lOKO Ce UZ = co - JZZ 
DDA „a O 75170. »»xX mu < ιν «αι 
ZO „ade >< >> Q >< Ni ᾱἵ - «О ч Ô > M>W Ω 
(C) — < x WWF WIN AO IO at ιτ 
Zw POO UW ZZ »-ZI TO m Z I UJ © €T ὅ οι zZ 
~ O сот. <{<dediebdeo Oee NO 2 > CZ em 
+$ e>OZ at ~ = το cq αι o» LU —Ju OZx==— 
- bk шж an zz zz zZ zZ, zou 0 tL de = C CS 
A In EF» — OOOOOOO-uQgrÓrÓr7d«u nl c - zuuo< 
at у>. ии N = x zc X — «t cC UJ cx 
ww ww DOUT ~w ccc axo ekea о < oc bk < >» 
Le Oe Se uuowuuuue WOOF - CC ucc aA vo „au 
„н ш we JSS ee JU TX: UJ << 00 DZ 
euj ә 53 —-2>< ы. Ји О Ја «то X OOo τ: ιαπ 
zz о о O ОО Оон Оччо OM LJ ш Ели - ο ος) 
uJz ooo ο e or a a EEE ZW Б Pr" LLUQ < 
Nerı an «τ ш WU te m Or uu ОШ ZF Nw uy 
— с) med Y <{ «1 чачды. чаш шм) eu — με» Le eT ~AT 
L N WU Tie tol ο YN O σι. = 2 
pi rd NOA Z° Lj XQ r X KL ТУСУ gd — X -—][] <f 
Z Z с “2: 2300N0 Aea DNZ > oZz „з UY) Le, c 
шшш PSU SLL UL UL LL et LL. SO сш Witty LI WW IA (^ 
Eom U! bee J Ωω). σσ Ξι- ο uz ue o 
Muy Wi Ζ O. οὐ κο” nn ο м хо = шола ο ιν πω νι & 
NA JA > б< СШ 4. JO μι EDAX ao 
O 10 UU) xy e Zum See zz mie © > XO 5 
Ст СС EU c LL au ccu c) yo) PI 
«T «t O IO = te we ο σας Z 
ο)» > QO ZM>0 G.G O жо) ке m 
» DO ШО 7x» => Cert | I t «t 
-- eC „> UY) =. 
= cN O ORI O MIT OnO Z tj a 
UJ e Atn O O eN mC) ese «1 wl O 
ΙΙ! 1 | = Ji | LIII T 0 u 
и Оч m DON 200 10 x m 
с? Ч ED OZV πιο uU) а. 
μι ceo o <f < 


QOQOOOOGOC(COOOOOOCCOOCQ(QOOO]CUOOOOOOU(OOOOOOUOOOQOOQOUOUOOCOOOOO(Q(QO0 


57 


CANOO29O 


IMPLICIT REAL*8(A-H;0-Z) 








ооо 
MAI Mm 
mo y 
ооо 
ооо 
Z= Z 
і< < 
CU 


ooo 


CANOO3O 
114) ,NBC (217,2) ,NSTRESCANOO31 
CANOO32 


(32) 95ELCON(105,5 4)  TITLE(CI1O) 


N 

со 

- 

E 

= 

O 
ZU) 
OZ 
cO e 
«Ὁ 

eco zz 
e O, ral sod ctm ame, 
Zr m» 
(e) eZ "uu ped — 
Сс e ουν e 
wel «---νγι 
OOM OW~ 
< оон 
QOQozaeodo 


AZ eZ wot 
L eO LJ 
eeN A 
~q DZ 2 e 
+> Z ФМ 
MT oe td ed 
ο ο] ο Ὁ ο 
`+ у) а 4 
(0 7» -— ON 
ST uf 
ч = an O 
S j= Luv cJ 
Cu) Z= O< w 
ОО 77 эмм 

wm, oo SI UL) 
abe AO, C) 
Qe A DAZ< 


NNI NNE JN 


ооо 

Suo 

anno 

OGG 

ооо 

ec s 

«(<< 

WWW 

= 

ÉS) ~ 

o + > 
UN - ---.. 
ee AC ΠΟ τω 
Ln = “UA. Q 
О = Z Z —=— = 
cN = а = = |j "= 
Ον aO je ~ sd ei T 
NO mli Die Zi I Du- H 
eC Па 7 Онч MOC A 
POO — Z Z O — <= 


2212220) < 0 ιο η ο 


OQ -*DIZ > AN =p* »>"SODOW NIDO=20 + OQ — Q O un GO tO 
ESOS ZE > JO NO OJOU 


=] 


© 
= 


ο. ο «οοσοοουωοωυς) οοο ώωώοώοοοοοοςο ο ο 
СО ОСО ч с\з Р (r OP DRONA 00) d UY ο ο OI 00) FIN OM OOO 
Q. QO. O, O O O, O, O O O. O, —=— e re pd ei e t e N (N (NCN CN N ONO) 
O O M ч pued pt — — —— — m ad MÀ —t — ed ped ed ered ed pad nd paad — — s s= 
CO OOO O OO COCO OOO C CO OOO OOO C O OCO OMC «Ό ος ο οσο Οσο ο 
E A E E E E E E A E OO OO ALOE ELE. 
«[ «I «I «[«I«I«( «f «f «c «c «t «IL «1 «f «I «( «f «I «f «1 «I «1 «I «If «f «Tf «1 «f «1 «t «1 «t 
DODDODDOIODOOOSOODDOIOO OC OOGDODOODOODOO OO OODOO OO 


e 
U^ 
U^ 
© 
~ = 
— 
-- О 
onl со 
> «-- --- 
2J u 
ι΄ Om 
E @ we 
Y O οσο 
— O ша 
= e e O 
= ο eui. 
w DO < 
YN O nO + 
co E= — о «ть М е 
<1 Q' c — о 
= C и 2 N re. 
T Coe) Б С O =” w 
(UZ = Dol w (5 5 
al roc) AX WA ο gr e 
«ca πω à) AYIA шо -—-—uC 
eZ O «eZOZQOZZUuUZeZZ Moe) 
M ο ο” ο υ---!! -- м оС > 
90 JONN eN Owe nae OU I 
"X C et UMN UCR he ee I QU) .- 
wed LT Det y e LO meer Yu OD 0 нч 


aw Se 


C.EQ.4). C4. (KODE.EQ.6)) GO TO 600 


і CO Da et 


w — + 27. C II ll O —G II Z Ji << O —%@ =Z n — O à = OO CN Z C53O — Ф О 
со UujJO T n Ox Ox 
NEM Uu up Do e πι ο αι m 

ZU 
QOOcOO-CuJ722zdgl í m OC < tU Q. OCO =o. O x Om OmU. Y JOU «t1 C07» ОНОШ. 
ош “оо Сс) Z Z £ rO. Z OZ CZ Z mO Z= >G Z CO Z — O <[ Q Ώρος) Z C— LS wt ee et 


© 
N 


= LU) 


O 
© 
— 


58 


ο τ ον 
GO 
Nm 


O Xei O 


о 
e 
try 


II O — — O 1 = 





OOOOOOOoOoOOOOooooooooooouoseiuouwouoo O O 
=A ++ LY O A O QA O — (Jr ++ uY ο Οσο MENT οτε ων e mum 
CO) 69 (0 0 00 00 00 en O AF AF f... r FA FAAK AAA AAAA AAN Oo № 
іф at Ч 4 Ч а qu med pur pal rl pl pro] rl ea, rra quad rl. med emt ed end ged pened nd end ed ped gol, con: pura d ч „ч 
CO CO O O O CO O CO O O O O OO CO O OO O O OO O OO O O ος ο ο) ο ο 
RZ ZELL ZZLZZLLR ZELL LZESZAD ZZ ZI ZZ Z Z 
L e LL ited L <T <T <T <( <í <T AL A ALE ALE AAT e < < «I 
00000000 00000000D000000000500000000 O uU 
—— aprum, nmm, 
м мм 
ES к е 
= ο ον 
> >- се 
aoc 
ο = > s 
Ф ~ & $ 
x >< >< = 
N NN -- 
= = @ 
e 9. м = 
= = o 
> XN H 
c 
e 
= т э 2 
es e p < 
-- x A a 
4 e Nett = 
9. es 9. Φ' a 
= e — Ф 9. 
m ο ~= — pam C3 
O = t zZ? ~ uJ 
O O IDO E uJ 
aul ә ™ LL) > ~) O. ш 
ο) "= ος => ы 
d» i mi ο ow œ = x 
«i «tr >=) ” ==» Y O m -— uJ 
az = ὃς ——» М Б 
---- O ο) ec cn ~ O o 
oo ONY < > = Ae» < - 
- + + oO М —M ш e 
- С. ч <I s == Ὁ oe = 
e =ч e - C <= = Ulm = 
< Oo ы = YW WAN T e e 
о б ~ Су æ eA су O e o. - 
— eo О <1 — Zen “OO ο N Qu 
«x 9 < O = ш ——— 2. 7 w 2 
+ © = O =O see >) = oc «oO 
= ο 1 < — «T LJ = COW = ~ DZ 
N ei © | C Oty << Y e N Za 
e Ὁ- + < =. - “0 шш © O 2 eco 
— =ч «tu» < 0а БЕ A ( 
~ QA nO c e eO. co ΠΤ) (ο 
O s JU =. e (DN (ut мч Y U => 
I mma s *N O = CN Cg co =) Ma, ~I 
O иц ἥἃἎ«“ SW COO α m e. e) s= O_o (Y "Az: em 
d el. 2" AT ιο" OZ ө Дел) ο ο τω Ow ος. O 
O мач di Zeel Or Se» oO N Due D ANN Oe 
1] u | ЈО — Oaa O MODO € eOe- A AN ο 
e-—C СО) um DIN IL N N >< — V CA HS DA e A PA o 
a O С) AO DD Al coU) eu) ro AAA ο) LIU οι ce eO 9S 48 ο 1 СУ mu 
ч ОО) + <( < — HO i i | + Doo e YOM OQ ο e Ü — — — O 


OO ll JJODO>— AMN UO dt —u «t «stueuoguaeujst«tu«zusd Owe 


«Id «xg <T < << << — —.D -QO O ——=t— + II — SE ZR Ee OFS ER Z= 1 = = —— = = o. 
Y Ur «αμα μια.) о О Оса 02 


O цан поо «i. p 


CO O οσο ο 
OM SG CAO ч 
Q O O Ol po 
omi omi mad pond ead — 
оооооо 
uL С 
LLALA A 
Saas 


CANO1720 


B 
N 


AD: 
YN 
OADC217$32) 935ELCON(10, 4) ,TITLE(CI1O) 


E D ME C: NCON( 1825 14) 4 NBC(2175, 2) 3 NSTRE 
, 


CL 
NT 
MBD 
‚CL 


LL Cx. f CJ eC 
Z 2D) 
— — — == δι.» 
= N U Uu. 
AO uc 
OO »00 
ora x 


-d O J O M NS —) —J O < O O LU << K O —— Z Q O LL £ O z O K O O K O > Cr O < O C OF Z 2 == OZ QO 
<( (5 <( (Oi < < << X (O(OC с K O Z y — — цц чо μμ πα Z Z u шо Ец ош. "OL 


© e © ооо ο κ ο ο 
`Q © N DD nnoon ϱὉ 
in © Ὃ ROO οσο OC 

Qe -- NA 


59 


e 
αν 
un 
ο 


OO 
OC 
OL 
00 


«1 





ӘӘ ОО ОООО OO OOOO O 
Ὠώιω ΟΝ-. Οκ NWO ὤΟ Ον us 
Pro Po Pa mm (9000000 DO DOAN 
en 
COCO COCO COCO COCO ο 99955 99 959 
Ze ze ZA (ZZ RE 
ΙΙΙ {ΙΑ чічч<ч <<< decide < 
ωώωωωώωώως ώς VV ωώώώωώωώώσω 
I <= 
eo 
2. Ele 
c e 
= =O 
z Ошма 
e “е e 
> IN a 
ca LL OC mi = 
= Or ο 
Z Z= O 
s COLL uy 
e UDNIN 
μα mz => 
— = OO 
Z ο. δώ 
e. OZ me 
α. a lez 
(7 «JO On 
2. na co 
9. Uh ADO 
WV) с Ошо 
ш =i O WY 
c (> ou 
= Ze DD 
WY к = Z <í <í 
Z UN OO 
9» σης - Z 
ο Q. ~ μυ. 
ca ZA OO 
a el CO) 
z Owe ο 
б» Ie (NU. 
O OM Rm 
«t a) ems 
Ο Ur» 27 
=j =Z I! = 
ы ° ο) 
z O PEZ Je 
e С ОСО тш = 
а. A — O O O + O Z 3 —— 
O ος <oon Ζει ον. 
к= ш СООЗ +. UJ” = 
и = чч > Pie n 
N ο) ο. OQ OD» 
- o = Η.-ΟΩΟι-- Όι- Σου 
LJ > ш < eh (9 ec = 
m е OC zt O FO JUPE. CX. uH 
or UO = У * eee S37) 
С ° - O ODODO e ZEZE Mu 
OA y s por zu = uJ m mm DLF 
O m Lis 9 SECOUS wu Zune em DL De 
Р-и е ар о ~ C “deZ DILO A 
N LTDa ea NN οφ М} „> O D шм)" ОШ > c= = Z 
- C C2 «X M MA || [} eee CODO m Μι) ο D oN Ze 
М СОФИ ПОЕНА но Ὁ mO 1 ı «e Z οἱ. 
D UVa rm OS DD О DON an SD 
YO m e < J — Фф s Pb "> e A DD <í 0 Qee . 
το γω j= CO OO >— — CO CO — OOO ο e ~w e ç ob (Jos € “Lt 
O) eren fe det LI LT rs e O O) 5 hm Pw et ee WY æ= UN md πι. O pu A LA TO 
Oe <I--OO-OO0OOOZOO LU {w UDEO A U N e d 
qi AOZ ww —— < A NO Ze O E O2222Z2ZPZ ο ο 
NE AO IL OOO JDOO JOOOX w OUulJUu.u.u-uu.dgo o» i 
ANA XUL -=-DODZODAODUDOZZ x Lime ZI UN LS 
(yA 
O O © O O O 
© u^ O ~ O C 
© © © O “ч N 
= == — == = s 
ыы. 


60 


READ COURDINATES OF JOINTS AND CONCENTRATED LOADS 


oOOOOOOOOSLNOFAMNOOOOOOOO 
RON A SUV OR ООСО ОУ ОУ ОУ СЫ С dr Uu «On co OS 
C O O O O O O OO O O O ——=— e e m e e 
NNN NNN NNNNA N СМ СМ СМ СМ см СМ AIA 
CeOOOOOOOOOOoOooooooooo 
E A AA NA A A Е 
<( <( <í <( <í <T <( <( Z <( — Z <í <( <( <( <í <( <í <í 


, COARD(IJT, 1) )JCOARD(IJT, 2) ,CLOAD( IJT,1), 


O 

© Су 

uy © 

N © 

“ч oo 

Qoo 

C б 4 
= — < 
= ο 
C Ga: C 
С LM C2 
=O ~ ἐς 
~ RAZON. = 
Ceu — Am ( = — — — 
ο -- » Fu < 
OO СО em Jin 
warn ' (rN — O (NI C) — —T <í 
еч С) Ὁ) ll O = om м нч 


Uo) eO cU mM OF OS OY) 7 — 
LUI OAL =< — 1 Y — 
Au O CO w әш) ү 
= „лоо ООСС < 
ὠμιοιμ)ή = AZII 
Z — i— SO TY ff aaO 
—-—— нч ТОМ И 
Uoc cc OU 0-10 
= ZOTO OQ -€10 Ux 


v= 
СУ СО 
(NO 
NN 
s o 





READ BOUNDARY CONDITION CODES 


Λι -ᾱ-ᾱ---νσςοσοσοσοςσσσσοσοσσοςοσσοσσοσσσσσσσσσσσσἒα 


CUNDARY CONDITION CODE*',//3535X,* JNT NO* 6X, CODE'// 


era 
E CN > 
N es 
9. e 
- H 
MIA 
-” e 
m 
- 75 
= e. 
ο. má 
pt = 
— ο 
O со 
a Zao = = 
CI ες s 
с Ом = + nm 
-— C a 4 mA C 
Orn Zea Oe. NZZ 


O € =-O == Olla = = 
ON O UN ret rd] at 


с^. l| Cw Y =< ἰὶ 
eM COO RAMO Y ч 
Оч mr DUNA, 


t Ou κκ — (> СУ 
be O) — po Γι ιο 
αι {ο «Luz mOO 
=== cN GO >= Б 50200) t 
mE Axm 
пс) Сс оаа © 


Uu ru Z= 
e SOS 

O CC C» 

со Об 

N MANN 


Ooo —w0 
OCMPLES “uLy-.-OP— QO C: O = CN О М со СО а ОМ со G: O = (NJ cO UN 9 Pa 00 СО πο (Y и 
смес сысысы CNOJOJOJOJ C009 0600 00 00 00 e eO o6 en «d wr S E SF SRF SSI UN LA LALALA IA LAMA ооо оо QO 
(NCNCGNICNCNCN (N (NCN CN CN CN CNNCNCN CN CN CN CNCNICN €N (N CN CN CN CNNCNI CN CN CN CN CN CN (N C(NCN (N (NCN CN CN CNICN CN CNJCN 
ao 5ςσςεςστσςςςσσσσσσοσσοσσοσσοσσοσσσσσσσοσσσσἒσσἒοσἒσἒο 
ENDO Mew Erg LER ZZ EL LE LE EEE RER EL LEE DEE ER FELL LE 
«[«[«1-12.2 5. «[«1-«I1«1«1«q[«1«1«f«1«1 «[«[«1 «1 «(1 «1 «1 «1 «X «t1«T «T «1 «1 «[ «1 «1 «1 «1«1 «T «1 «1 «1 «T«f «1-1 «p «1 
ae O LE eet Ce OO OCONEE reese 
τος 
~~ Q. > 
@ ο. P= 
= a 
ш ш 9 E 
е кє O = 
< < =ч © са 
2 Z N = O 
Dae} fond > 2 
O Q E Y Е: 
с с o ο x 
Ὁ Om (9 
De IL = = 
O O Oe э e 
at © x 
> c O and > ο 
E on ec» LA 
δ. Фф» -- Φε 9 > 
< ὃς 4 >< O > 
n ın m Ζ O 
& e bad ο. 9 2 
ш ш O e <Í --- 
= |-- <-> >м ш 
<< Ο- мч 
z Z -JC5 = ^ 
= мч O = > 
0 O ^» === xu = 
c ο UN + mu ~ ~ 
O O N e e eu) М. + 
о О eu e ч“ ч e 
Q oO pal ө Н У м = ©» = 
e $ =) ος) oy = «Y T 
x N Q - em dl x "I = 
= a ecco — г e =ч = ». 
e ο. <I » =) U þe e ο ο. ui 
<= >< Ox ο. rt i = u II 
in ow Om = = = dq. 7» " 
e ee ο. & <I ce - = es Bs 
æ 5. om | СЭ = LLJ E = <= 
c c vt 0 pat a ==у > > = 
UJ 113 am aj — O e = к= = 
СО — —— Y) Z c — > — ~ -- 
ENIN —(5 ш e a 0 € > = -— 
De De O — O <I Zo (С -- O = 
2020-6 XX br =) „I OM Б = e © 
< ἴ «tu чо c ul σα Qe < ce Z ο 
2 СОН Om ш <= wd цою zx ш == = 
ET у = Q eN XI Wt 2 ο. ~ a 
v) mM rn - mn, m k= e fp nn Ma т tL es LO 
Е бои. ы X PIO F κ ex + O “> ae 
m "η η (С = O. <> — < 2 Lu кз O Se J $ 
—ч - в Hm x με 2 “e > -- uL UL) m — 
^ a ma &FOZOW a Z TO Zen - Ou amen =+ 
== e» eX eO «A mn ο ο ο). F Oe eC HOM 
η МОМ С ri чч ΓΗ lee ND O OAN et Ot tf на 
=ч “Э м ο  -- Се l Oc. I ii cN = U cN < NULA HOY ο» 
a? Neo Ne en Bb MM (OIN AAN С = as, pt (N) —4— e >< 
Due Due Dom O I— — веб ау 0, v u^) > м) а @= w O CN 
їг ΟΡ Όι- ο e Z омы н он © Б оно LE 
O ma c L) <L m) T OUJET CO ~e <L UJ <T Z — UUJ μι {πω «ουμι«{ 
Oz ox EZ“ С ADIPZONMSZ CO кє CN OC Z CN i = 
Н Ζασαρα μις <€ ICN ια l wN ANA € 
ο ο  Ὁς Ου ω ο Шш OUWOXONOXO ш αΩο Ου οοασ ο 
DULL aL SFOs c OCW SL συ A SUL WOW LOL. 
=i = = --- 
GC C Do COGO κο ec) O τ COGO 
Om Y IO Oor O сэ сэ © DO. OG 
+M M UNV Oh DON © "oy ea Su or 
οί - ο αν rr (N Now NI сус) Na 
OOQ 26020) 


61 





leleXleleleleleleleleleleleTelololeoleofele о о e © QoOoOooOoO0oOoOoOoouo0oo 
OM 00 Q^ O e 0J rn + О ЊО О О —CJC00 «T UA USB € Ww > OQ —CNIO wr uso 0009s O — cN r< 
OO OQ rh MR r= Re EP = СО со со со со О со o Ον Ον e O O O O O O O O O mie ed e 
ANN CJOJOJCOJ CJ CJ CN ICNCNICNCNCAJ N €N CNJCN CN CN cN N N N с сус СС MO r MO MO MO OO eH 
0000000000 0000 OOOOO о © о © Ооооооооооооооо 
ΕΕ гг г г гг Z == = = = = Z= Z Z= = == == =G C 
«T«q[ «f «p«r «I «c «f «€ «r «X «1 «X «c «( «f «f «T «v «t «T < < «t1 <í <( <í <( <[ <( < <( < <( = <í <r > 
QOQOUO)OOOOOQOCO OO OOOOOOOOUOO ο ο“ ο) ο Ze 
— К 
+ T ” uJ 
-- + e .- ш = 
$^ >= = ri Cc 2 
= gum] qa Φο з = 
II e II - >< = 
- » E Н ш — — 
ш II «t ο O 
= = e è ~ і 
= = < O MS - 
> — О = - Uj 
O = - - O N - 
ш e ο Q ⁄ | 
e = ш ze OU = 
ш O a 9 e k 
ш ш ш - O c e 
ο) e uJ 9 2 = - -- 
x LLJ (2) O = + ~ ++ 
ш ш x x O Q + οσο) 
с) uJ ш eo » ή Oe 
o >< + ~ — 
= ш =) c) ы м w N м. 
= <ї c = no o >» 
ш VY) 4 < t= Ф xt OZ 
= = α΄ C = nn e OW 
uy = ш - + > ы” 
--| рае — Un N O im 
=. ш O <T 9 = ος) mf 
< =j = Ω Ζ ΤΟ N 
a LL, M NO *. CN e 
¡de) O LL, LL. LL, O - ας) e CO 
O О O © =) Y w* Z m~ 
r. ES «€ MM e w 
ο Z m - =a a us е 27 УО) 
a e = = - O mn) ο” -- 6) 
ul 7 e o. es - +O 90e 
Q E <> = e ο e = СЧ <t > < — 
о -" ш ZB -- = = -- = *QOO+ 
© > 0. QOO e e =ч e 9. TIE. = 
ο ος + + = > e »” e am QNOOUZO-ft 
e ш етть A tus = am a rd) Q — = >= 
NO > AN — ”- Que = Wm Yt «О m He 
© ° J ^ ^ α΄ < e ZN eb Ju) 
= = O Omp ul e Me OF Ca e | Nil eo 
о ш ο... «ας к IO JO Q~ „JO Zh. 
O qc ο τήν» LS Cow δα OG We SZ will 
[9 — Ы. = 22 эмче то πο στο ο πα Z e Ci ee» e 
[N Z San “7 az 27 == eZ Z= | + — N — = 
= Ze “=ч „с а= 20) = =. — weet AMANN 
~ Om 4 X e *000 ule ο Se Me De Цу FL ου 
9 M+ + Ze aa) J ο ο. > A ICI Уус Lo CN Το) 
Сї Q. > — к= + —)w— =, у) < <ї <1 < < C st vu ne 
ш Zo РАЙ мы O >Z Uli — ME eu eu ο Wed Ou JI> Ὃ 
° "DD N II М (OC) ee——— O O O O O «f =< m= = + 
- mM) D m > e C —— cJ Ne fie Ne Ne Ne »- ΠΠ Qe ie - 
ЗО 1I Z7 77 — ( Z IIO = == Me te Mee Oe & am ш: LIke CU 0 — 
O =H U » a DO + ἰ--- э ч dese = з» Y πο τι ο ο η 
Dr aw (OL < "8 in min mia elf et^. 4 — u) U) — — LL. C) — C =“ 
C) OO ο ο  - Pole Zu < E c bm DOUNZ Z JZ ZÜ eO 
ODIA”? ZZ ZU e«<ADO-XiIU UL Wet Wet Welw Овч 2 OC, =OC Z OO 
Zi ZO error Y ++ Z< TO Ба. но. oO. = MOO α JUD = O. DA lI — 
~ Н VOUDOZ Il ~ «= сс) —gcO0O-oO0-—o0oO0--gc Онна Оооо. = 56005550) 
LOKOMMOMO ш шошо Јошион- хон шон шою хог 21020 5 -оОО2 оошсо 
= бу ЖУ) с -HLULOOOTZEZELNWZLDELNIZUNDZEUL.  MuUN-=000 —"2)OXNO 
= Eel 
© O O = Oreo” r 00.00 
O © © © O TOA ο ο ο α ο 
© eo © = nam vr NM OQ rn 
vy (o T + + Par tF? ES [e FS 
OO (C 


62 





ο @ ele ok iene lee le len ml @ (al ela ieee elokelele la |slelelelelele lela ele! slelele) ele) ae seal ele 
FIN DOM WHO SAM FIVOMH GO GX O — (NJ MOSM OM DRO SUM ENDO O C O — cN r UN O Њо с O ч 
Mrd e — — — (N (N СО С СУС OAT A CUM CO 00 00 00) 00 со сосе PS RO ORAR LALA LA LA LA UA LA LALALA O O 
co c A COMO COCO COMO mn 00 00 06 00 COCO 00 00) 00 00 c c 00 00) 00 00 08 06 00 00 00 00 00 00 en 00 n Q qm 00 e en e e e e 
ου OoOOOOOGO]OGOOOOOOGOOOOooooOoOOooooooooooooooooooooooo 
Επ ει ΕΕ κα σα σα ut ac == =Z = Z= == =Z xU M ασ ο σώος σσ σα zZ 
<í <í <( <( Z >` =" Z <( <( <( <( <( <( <( <K <( Z «r«[«c «t «c «C «C «rc «I «xt «c «c «E «T «x «E «x «T «X «x «c «c «x «t «C «I «X «t€ «1 «t 
QOOQO(OQOCOCOOOOGOOCGCOCOCOG(OCCOCOOCCOCOGO (COOCOCOCOCOOCCOCOCCCOOOOCOCOCCOOCO O(OOOOGO 


e 
I! 
= 
= 
k. 
v» 
v 
LU 
Y 
= 
IW 
= 
<I Ww 
ESL. 
а. си 
κ.» 
σι 
Oz 
u. = 
> 
xo 
— με. 
cC x 
< 
< 
=v 
© 
O 
--- E че, 
v») A EO 
uJ = WY || 
or ~ «422 
шо? «αὖ LL me 
cz =] 
Be e wo 
ma. i x τω 
= < SZ - be LU 
eZ ~ = Y 
== Oe х. zov 
N X Οι” O = mo Ya (Y) 
+ m Or» O N cr LJ 
-- e oco co + а. і ош 
- N «5.22 e O — ο)» > 
{ — ri a COPI U "We 
— o == I zxcO Н N η eo .J 
me C3 QN ου ο e arma 
ES Y © h += F (> w ~ f <T 
мм W D-RO mhna A Zu jJ -- ~ — O WMOWO + 
ulcer © 09436 Т нч О) ОШИ + ` σσ ο «3t 2 » 
+ Ζω Z NII ZZZ eJm * 000 =D C5 το 4 
N e »= 27 e ALM HUU OZ ZN YN" In ZO UUJ eu 
O — — -- а ted [du u ο ο > = + é | ο Н = % (J ΕΟ + 
-CCoc OW Nail Or SOL «ΖιΏ Omen mm Omen rt (D 4 τετ. - 
BE men ο ο το ο etit tcr» woos = Pez = z ed 
SAW NE meme 1 dm Zu LI TU KIO NL 1 LÀ ш Ξ me C CC O | 
~ па. -- mOwOWOY d: atu --»-- ë Dn DDO => pæ р) в Гука 
O73 00D BR ZO O~ 220 NOR WOwr C00 Om HOS Pee LODNZ Dy) оо 
OO ο ολη) 06J7 0 ΙΙ Ə 7 О 30 WY OOOH OZ MIE RS HNINOW I Le Om Ze eO 
pet et ONS rene Sae I Y фр 4 аа СК ZO) IA LUZ —J —— — (N CN CN JU CN Nee pe HHO 72 KC JO a U 4 CN CN 
м.м кеч I! н Ок iur Io _J— e uo I“ I LZIZZ ZO FONA ai 1 


Qu zodoco»zxugogcwocco:sti zoxu.«t «xou. 0o] -3]7n^757u JOOOIOOIOMMUMUZIZTAHAZZO 
QUUD OO J JZ C — Q) CO TZ ROAR HOO Ewe ODOM YD MH CcOCc)0 09 OU «t OC uu P UY)XÀYQO F- O 


O © O ° G 
O — e ° <) 
== = —i e e 


6З 





@ O O OOO O O O O O O O O O O OO O O —CN rn ++ 5 .OrS— Q a O = (N rñn+YyO O O O O O O O O O O OO O uN 
cN r + ОР со СУС — Cu ++ Ur .O — O O° O O O O. O O O @ O O —=— — — — + UA n Οσο πο ОКО Ro 
DO DOOM xQ .D r = t— mm m me m 09) LO 09 00 c со СО со со со СО СО СО) со су со Сос) СО O O O: O G) O SONA ON O. O. 
PCY MAM n n n nes nc 0 00 e 00 rn e 00 en r eye rn rn 0 en 0000 0n rn 00 e en e en 0 e e e 0 e n n nen n 
с O CO OCO Co CoO CO COCO COCO CO CO CD COCCO C OOOD C оосо осо ос Оо 9 9 OD 9999999999999 
ΜΙ ευ = > = >= = < г == ο σος CCC οσο 
<( <( > <(<( <( <( <( <( <( <( <í <( <( <( <( <( < < > = >` >° = = > >= >= >° = IT << <( <( <í <( <( <( <( = <( <( <í <T( < > 
ο το τς το OOO το 


Un» 
Mm 
LLJ ame 
cc — 
ho — 
N) um a 
— Q. 
zZ N ως | 
— YN OQ. Z 
< LLJ CC +£ ο 
ο a HO u 
i= = шс) + 
un VY) =O © 
ze c 
ш uj > см O 
27 = nit á е 
< < it N 
- a > O эшо 
= + — 
Га е 4 ο ος--- ος 
Qa. O DE ~m а. 
LL. tt LL. LO m х 
= = + Cov СЭ 
<b --- XK e xZ On. С) 
| wu sx HO Γι © 
E tum ~ d mn  ( E 9 
FO Ó = C me 00 e 
eS cf A, O <— e. < ww | 
e | S N i κ s = 
і — Z очы еа À 2t -— — 6 — СС си — — — (η 
OO WOLLT è OF αμ (m Ссс οσο w 
MO ——— А. = тым er Le mu x= (í. — 5 nn we 
һо Ф ——jj$ Oo He x«l ας 45 - “νι κ 
- α NZX «ος. ALA NZX ee Инә ү Сү 
IA Uy ϱ το ιο ο et 2 uuculcoco-aouiuJ 
-) i ——J-— 4 n»n Hos d ән но тон ои aS ou HO Hoi o Hu d 


ONZ νο ποσο NNSA WHAM 0900 C400 «f cc 
СО eS mm tr 0 9 ODT ο ο ο ο er Dur a a mmm. prom) 
— — CÉ || — — (N (NQ) = r. И ч СУ Ју ОК ч AC OJCOJ 00 0n en Б О 


CO v-u.uoco»ogocdoaiuudocoococozocdcu.uooooooqoQqoocdccmu 


e O 
e = 


c 


e 
m~d 


e 
= 


64 


N 
N 
uj e 
c» 
ц. 
ц. ~ 
— = 
bct 
ALL © 
LL 
WIN 
Τ»-- 
HI 
bh UN 
л =ч 
zz 
Qc rU) 
O U 
u nr 
he 4 
ke Z Z 
<L LLJ mug 
ILZO 
f= LU сс 
JO 
LU LU 
a UJ 
η.) > 
< 
ЭС. 
ου 
a ru. 
oag 
AA 
Ne 2 
СЕ) 
FO =. 
O 
rt >< 
ime CP af 
«т = 
ev WwW 
COU eT 
WLI ar 
Y TIT 


N 


THE STIFFNE 
I 
E 


NTAIN 
L 
D 


N THE FOLLOWING COMMON STATEMENT 
CEC 
ROGR 
US N 
ENSI 
ISYMMETRIC PROBLEMS) 


HE FOLLOWING STATEMENT: 


i 
GP) 
I 


24 


O whee SNS a 
ujziuuom** o«u 
JND ш мч 
-= SZ J> Z Imee 
«Ὁ Εμ) ΤΟ. {α-ω 
ONS eqo 22 
O kh Ls a LL) cem у? = 
л м сс van 
PP ICO ο» MeL 
TINTE NN 
LU xa QC e «1 UC zT Tr 
ZO Бола 
e d) C) c < Y 
EUNT) = 
GANG 
ο σσ κο L 
COD. = 
> ¿< Z= л 
n UO 
c m 
ш ш> 
жш ү Tr GD 
GIO 


SS AÍ 


QDCOCOCOCOCOCOCOCOCO COCO COCO CO 





C O OO O O O O OO O O UA O OG O O OO O O O O O G O O O OO O O O O CO O O O OC CO CO OOO O OOO OO oO 
ODO SAIC «f LAO Pa Pe сосе ОМ O Gs O —C(M rn σσ — “(СУ „е СА) n ОУ O MANOA 1С О К> СО ОУ O — (NJ 
OQ NOOOOOOOOO O mM M m cm — (N (N (N (NJ(N (N (N (N CNJ ur 0) 09 CO OO e em e — cO (0 00 0 c en en «r dT + 
VO a e Se ee re С rE rr rr rr Er EMO ES rr rrr er err rer 
@ O O O O O O O @ O O O OOO O O O O O O O O G O O O O O OO O O O O O O C O O O CO C OO CO CO O O 
2 22 721222 ZC ο ο ο ΣΤ τσ τα σσ σσ σσ ρω” 
«[«[«[«L«[«f «1 X <( >` =° <( <( << —<[ <[ <í <( <[ <( <( <í <( <í <( <í <[ < —<[ << <( <( <í <T <í( <( <T = <( <( <( <T <( <( <K < < < 
ο σος 2. ος OOOOCODOOOUOXUOUXDOOUXSDOO(ODOOUQUQOI(OXDDOOXOCXXODQDOOOUOI/UODOOOX]OXX)ODO0OUOUOXUOUQOG 
edi o 


OVO OQOOOXOOOOGOXOO(OOOOCOCOCOOOCOXOOCOCOOUCUUCOC 


IS THE NUN: BER OF GAUSS POINTS VED IN THE INTEGRATION 


COORD( 12,2) 


NGP 


2 
O = 
WNZ 
EE ul 
OZ Z 
JJ LO > 
=) а. < 
Ib NZ» 
AZOWM 
посо?) 
QOO πω 
O mad 
OM ONT = 
Шик шс WM 
Б СО си — 
< NEMO 
ZOIN e. 
— = y «ОС 
O wuah 
Ce >ш 
ozuer- 0% 
OWI a= 
VZI-UDN> 
ul Q WM 
ш ә сє 
тш Ох 
- «ul < 
ul A> 
со си >< СИ 
Z= >O 
— = <LL 
Zu Zu» 
JM CT =.) 
I «[I-»o 
EU) << 
2.2 «х 
OZON 
Del) ZI — 
О шос 


ELAST(4,4) 


ETA 


ETA 




















.-- 
N < € e 
= = E =j r= 
I ^ | 
y + { 
N + Q^ 
Em + > 
I * I 
P ΘΙ | 
с со 
= => 
t { 
I | 
+ A © ~ 

4 

>< 

+ 

+ 
e со po 
= = = 
I + I 
I + I 
I + I 
I + І 
N τὸν Ne) 
+ + + + + + + m 

I 

I I 
I I 
I 1 
a `+ u^ 
= > > 

-— 

x 
= + = 
I “= i 
і + і 
! + ἰ 
I + | 
I + j 
dd 
| I 
і 1 
l I 
I І 
N 12; 
= = 





65 


CODED BY PROFESSOR GILLES CANTIN AT THE NAVAL POSTGRADUATE SCHOOL 


~ 

Νο UY о 

Oo о O 

OO N ο 

Oe $ Ww 

nı m © 

Ц\ ~ . e 

no e" C 

nO O м 

~ мо O a 

un > к GO 

y исо N ο 

un UNA NO + 

< NAN о о 

e On + m 

-. охо ο ια 

$T altar «5 ° 

~ л cud ιο ι 

+ + HAD + ο 

<I СЧ чо „ч O 

e eO ++ = N O 

~ (00 NOM N © 

ӘӘ FO лое MN ° 

~ wN OD I e O 

«€ 20 Чоо” e ο 

< (CN οσο ο © 

e eco хос O O 

-- κ ADO MA N 

- - N FO FOr (Y LN 

Z Y амо Мосо O б 

ZN NN eN NOANA D m 

e ~ «< “OO eco $ O 

O. οὐ (οἱ Mm et 

c5 ~ т и СС) ο) zo eo 

Z ω ων Oym MW ο 

e ο er COMM - y 

Z OM ODD ~ %4 

A A е *О N со 

ο cC = = =+ | Oem W m 

2 Rar q φ» e e 1) О [Te 

e $e TNO OO • © ° 

co Nu u DD СОС ~ е. ” 

raw SOHO MOD o O 

MAIMEEIK™~ON омо ο O 

< | υ οσο UO w+ o 

SON ame ID int τ- ο 

L = BLAM KWON SIND οὉ ο, 

pr Tr a б нч Осии ~ ON 
Dior NOT MW WN 
wA NAGNAO ON s+ Oo 
New % e mm чм, лу} 2 сЕ 
CA OD f—— —. F Оч ош О + OO 
{ο ον ere DUNC uU S 
Ὃ.Η----... ee DANN сә ο 
Cr Ч-СЫ 9м. М Ον 
UJ«t»x»«»«aoct-Or-uNoO-* Own 
UJ cí MV eNO MM of 
2 ZZZTO ert ο OO eO 000 


ОӘО о SOS SU 
— —— —— (NICO UO er CA ο ο ανν 
> CO Q m O Z >x «x «xus px e 
O»=ZZZZO ec c y 
CC ООШ «f «T1 «r1 «t «T —<( — <í O; 
m o. = >: = > = = = ee Mein Q 
D E C) < <X << = -O < <í e 
Qy Сб С С CC CC CC CO OC О 1 

= ra = 





ΙΙΙ -«-««εσοοσσοσοοσοοσσσσοσσοσοσώσσοσσοσοσσσοοσοσοοοοςσοσσοσσσσσσσς 
CO АС ОР СО AO —CN 00 «fu Oh 00 ОУ © eOIe «tut 0 0 O eO ЛО КУ СО СУС н СМ! сс DMA cO OO 
SS + + ώ ΦΨιωιοιω ALA IA LA LAIA ANA 0 0 0 9 0 0 0 50 60 30 Pino IS f b Pw PISIS O5 0000 00 00 00 0000 00 :0 OS 
PF Pts Pet PPh Fh eS ae E SER SEA ESE SEE E SE ee ee tt ret rer ee eee 
коо ооооос ос осо оос оос ос оссо 95 95099509 99595 9995 959 995995099 99 9 9 8 9990990909 
ЕОС СЕСЕ С Ns e e E a A e A AA A AA EZ EA A e E 2 
LLL L L L L L <1 <1 <L <I <L <L <L <L <L <1 <L <L <L <L <L <L <L <L <L <L <I <I <L <L L <L <L <1 <L <1 <L <L <L <L <L <L <L <T <1 <T 
WOOO UO UO OO OOOO OOO OO VOY OO OV OOO OOO VU OOO OUOU OOO OV OU QO QO Q 


e ri rd } 
O | = = 
e с) О 
O OQO 
e QOO 
nei © fe 
O e e — c 
OO mej 
ΟΟ | О = 
о eco 
ч о OOO 


Im ο м 
=en Ommo 
OO . * nO 
an Oman 
оо rans 
ee eC) eo cà 
ast OD en 
Ι - Ο-π-Ηςηςση 
sO Cn =Y rn 
OO monn 
C) O < n Cy nr 
O sO roro ror 
* C» C3 0000 e0 en 
=“ РС) ССС) егеу 
FO eA? OO 
ОСО nen em 
OH "MMO e 
О ОГС) = 
ϱ- Ὀσ OOO 
= | O rn rn n O 
кр СССр сус 


78628670499300,.568888888888800,.478628 


OO" осууго © 

οι I cene ll С 

OO = =c cn y D, un 

e OO em) am ° 

ч С) у AO” - O 

e ONOMNMA 2 FR + 

O = en Cy nn а. x - 
wf OO N- = ~ — 
«ООС і сусу 2 x m 
eC) O FMMOMNMS e =з ~ < e 
Онч «ON MW > e = + 7$ 
СО і чесе amna ~ Z > -- Q 
AN =] ONNO e e u^ в t= Y (2 
oO» «MMO | с. (o z < = — 
тил OO CO — (9 en o e —J = = i Un τ 
Оооо і гсп OZ e = co <t + == 
O eO Mm De f c = <> = о 
ο ος a Pe υ ГЭ) Hd >- =ч uM СЭ 
Όσο | ΟΟΟ MOD ec =. <t * = e О 
CN € OD ec а С) e 3t >< - ч φ. > = 
СООО бо (= со n. о о-о. A num x x = α 
OMAOAI ea QO Om OOO O x5 — > — Zo 
Οσο «О eHOO2Z2Z0 > 22-2 Z < =ZZ"* we ei^ 
(V) 4 e eD eD rZ а m © » a gç m — — — ο — aran > mit f A 0 e m nU) br e o 
QJ e—— C) DO OA rm O а (NY td id ν-ν- нч на rm я =ч Mot Dre CA) н Ф 
ο "Ον Ο ον μη Ιον DO, Oo NN N TTS 
<WAOJO ος) ο) οὐ 7) —M ^ e ο СОР р М М ИИМИН ә на ee oem OO 
LAA > de te ιο >) <= > <í — > <í — >< <í — >< <í —” -ο — —C < -η Y in CN i= CD s 
IDK => | — ә ОО =Z O H I H || l JJ Il Z OO =O — O LL. JI OO "Ооо = = O. O 


Ον © > ο ΖΟ0- и Y+ O — — C)— — O — — Om RENN HOW Ow MOO `< O O 2>Z Z G O 
<í + <í < O) < O O C3>— = — < N N JO ALT er ||] me = 
= O i= — GO = GO O — wa меме a ee 2 <“ >x x 1 <t м мм O 
ΓΡ) {« οσο eu OO 70 0]0 ——O ——OC)-—-—CO)——OOQO*-O i O στ.) ΟΙ ο. «2020 
ООСС у о оч СО УС) Со 5X «105 «160 5« «1 (0x «1 0 O00 «1 C3»x O3» 0 «103 3) 03 UY oZ i DH 

= = mat AO 


° N ra ar LY O O O O 
o (η e m ү $ Uu n O 
-- ma m=i == =i «5 =} N eo 


66 





o)sle)lelele)elelelsls/eleleleleieleleleieiele'»]eleleleleleie 


000000000 ANM<SFLUID Mn 00 


HUM PUNO PM СО ОСО ч с) СС f UY O SOGO GN O = Ç(NJ r ОМ ОСО FR NINA «ο O QA O O O O OOO CO O 
O. C" Q. Q" Q.x O" Q O° CQ° C;S O O O O O O, O, O, O; — — — — Sh Ht HS HAO HAIN AI NN ANA MMM MMMM MQ 


«ΩΓ -Ρ.γ-«Ρ-Σ-γ-Σιωιωιωιωιώιώιώιώιώιώιώιώιώιώιώιώιώιώιώιώιώιὦ 
ао оо оо оо ооо о ооо ооо оо оо ооо еее 
ἔς 2 E ez zz 9995 8 2 2 22 ze 2 Le 2 EL 2222 
<( <( <( <í <T <( <í <T <T( = <T <( = 5-5-5 <1<1 <1<1<1<1<1 < <1<1<1<1<1<1<ї1 


US UVP US UA UV UY US UP LALA LA LA LALA LA LA 
00500000 COSO 
Сс суг г ccoccoooco 
<( = Z <( <( = <( <( <T << Z «12 >= = > = 


QOOQOQOO OOOOOOOOOOUOOOOOOUO OOQOUOQOOUCGOOOdOoOdooocoQgooooocdo 


M)tCLAST(L$;JJ) *BC JJ 4 M) 
? 
) 
T 
F 


KX | 1 3 B3N$, NN) 


— wm ο » ll sm et TO LL cx. « CX. 
“NO NODO NN O —<( ч І Hoo = -( 
= =O = = ne LIZ TIMO Xamit = he UL 
pu pd [М Б uo 0 e 7—7 ~ m TA). UL — — = 
u ww ον CO ZO ZC +O = ==Z DO O« 
J-JUOOODOOI OO «О υ CO ον Or с 
De >= = >= >= > >= > Ыс). с. Ον μον. скы Du 


xx XX XX) ~ ~| ee Ze OOO. = 
H HOI I OII I< OOO RKORnIONZOWZDI= YO 
>< >= (5 > > (5 >x > CO O0 TOT HOHONN UO LW Ne Z 
N< 
© © ΕΟ O СО = 
O ° ¿10 OD O DO T 
+ ut ο ~ co DO — > 
E 
QIOUUOOGC QC 


67 


21 30ONX( 23112) 3W1(2512)58(4,241 


4,4), 5SS(48,2%4) 3 SN( 485,524) 


= мо N 
«t «to +20) -- 
et) CO € (СС nec) °- 
Z—0 QC — LLC, OC C O O =] 
СОО ОИ е Ὁ | 
rent ya es Ц MOO . € o6; 7220 
N+ZWO æ о е о С О ММ 
20000 "NN ilz ч 
UO RH H CL — 
=. 5 JN Ос 020220) 
OUA UOO Z Ou. O 
Ο eQ-—nmOOgOOrrTuL:7Z-Z-COC 


© 
=з 


C INIT 





οο ποσώ ώο υ DO LA 0 Me 00 DO NNM ADON οοιο Јо ио Љо Оо мело 
DENN AN AN N ATI RA ARSS GÍA SS FS ωιωιωιωίιΩιΩω О 
PED 00 00 00 00 009 00 00 00 00 00 00 00 06 00 06 00 00 00 ro cQ 00 00 00 co 00 00 00 00 en c0 00 00 00 00 00 e 0 001 0 06 00 0 e een e 
LALA LA LALALA LA LA LA LALA LALA LA LA LALA LALA LA LA LA LA LA LA LA LA LALA LA LA LA LA LA LA LALALA LA LA LALA LALA LALALA LA 
DD DOSOODOSSOSOS SS 9988: 9356555955 9 985 9599959559 9 9 9: 9909588 © 
ου sa (zo =a ma ss Gal ss malas aja [aja @o ss am em GS Ge E Ge ms dm Gs um G qu G malla se s saam i am Go c sh “siasa dm as s = 
LI MI ο ο on nn nn > >>> _ 2 sn x= == «д 
= — 0 G O C ajo o ajajaja ajaja aaa jaja as ajaja 


< 
= 
LU 
e 
Z 
< 
fd 
>< 
pe 
c 
= 
Ф 
= 
O 
Pd 
> 
оо и. 
OO OOO 
O O0 OO OO осо е 
O ο OO OO OO a 
O NN + «е OO OO} a 
ee CCN NAGI 9 5 9 ο LES 
NN ММ AN NN NN == 
NO ame NN AM an = 
~ FH ++ << >> μι 
О ии шии dex ar τς 
c (Y ο тт TE CCX T > 
ο αι Ὁ B = me E II II 
Ge ЭГО + { + {+ mu be u. 
i © O u + wu wu + | + 1 ο 
V d ш cM ~ ZZ ZZ WW шш 
e— Мм Мм «м + OO OO ZZ ZZ ° 
ать om - - ο. ~w ~w OO бо > 
z zgouguozo = XxX >> KH HH ex 
O Ὁ σε τος х аа στὸ мм >> LU 
ο Сс Y = + OE 1 = i COLO == ZZ O 
e CP CCo cc = >*= | ` I <>“ >x ποῖα CO OO 
= HORD O + — > mm + — Qn — — > < Z >->- — 3 ++ — 3 £ LL 
©) Ра -- ο Z; rom м, х 1 Sch x ” <. ο μας ML Eo NII A~ М 7 
-- m | uuu > — “> Κο και ж CORSO no τα 
JD ο ММММ. и O» OO (0 w ZOH X 23£ 3e 3€ OCO OOO OO x 
С? + г —— am > ἃ Ge Fe de de de 00 2-00 Oo OK vv irn ma 
> er D >> >> ш XZ >> > > > > > VER In ZA a x 
LO OO аса ο σα Z *^»OOOOC»Q230.»-—223. = 
= n ο) cc NEL o σοι — Ne M36 4 aen. OO C0) «Т 
оч e LJ δ IL idee δὲ Lew CU Lu >< LU LL) Z $£ + C 31 3t = 
C) Жэ жь C Xxx A 122452 >< > >< >=> tt zz ZZ U 
OW tili iee d SAAS ew OZADLALAZTZrU U Ku Deu X ZZAZZIUI м 
ешш σΣ Οσο ἱ.). ο ωστοσο ιὮ Anz ZO = 
О ZZ Ζ” ων "uu < a tueri к= Ü] ` ww 2. <. Z =< O— ιο” σσ” ” [a 
1n 200000 nw icon m " l H H ll H gu "n Con Teil H H | o m 
=Z Il ji ij ir СС — —— i í ων πο] NN we RON ee em eee »ς 
— — >< >= >< >= « MIAHAN SF x uy Or Or O+ Z= Z 00g ——-———Glcosdruoz N 
Uv» ος Pad NVVN OUV << NIIT (0) UNI UI UI Y) UN) UN UN IU QQ) 
— = aa = 
emi sg τι Л © D m” c 
πι -- — — 0 — O 
u. 
QUU Qoo Qo WOW) 


68 





CO O O O O O O O O O O O O 2 O O O O O OCO O O O OCO O O O OC OOO OO O ө Дә» е ә ә эө], 
OOO — r Qn fun OS O O-nrm«tuvot-c 0 O e0dim«stuot- c0 O enm «uxor O o аме 
MENS SPSS LALA LALA LALA LA LALALA `Q QO O <-O O O <O.O Q `Q r= ID PP PP P = co co co eo cO OO 
UNUS US UU UO UO LU UY LU UU UU UO UU UO LU UO UU UU UO LU LALA LALA UA LA LA LALA LA LA LALA LA ALA LA LA LA 
OOOOOOOOOOOOOOoOoOoooooooocOoOoooooocooocoooooooooooooo 
22222 22 22 Z Z 2 Z Z ZZ Z Z Z 2 Z Z Z ZZ Z ZZ Z — =Z Z Z ZZ Z Z ZZ ZZ ZZ ZÉ ZZ ZZ ZZ Z= Z 
<£ <( <( <í <( <( <( <( <í <T <( <( <(<T <(T <í <( <( <( <( <í <( <( <( <í <( <( <( <( <í <í <( <í <( <( <í <í <í <í <( <[ <í <[ <( < < < < 
00000090000 0000 00000000 Y 00000000 OQ. C90 OO OC CO OQUUOCOOCO QU) 


a 
τ τε E CT т оу «т v 
O ου Ὁ- ΙΝ, 55 -Ὁ 
со Ec CO O OOO и жос) 00g 
=i ЛШ, АШ, πα UE үй гол әш LL tE 
= De nn A NEN хх ο NER IN 
= Е c »- FNS OI KKK OSX 
O ( > r > N O a Zz>a σας axr 
PZ e — > >x > = >x >< i= _) > C) > > > > > >< N 
e СО ujeuJ-UuJuJoe > e RODA te er rl 2 Z 
DO w Z Z Z Z< (C) >> >x > 3£ >= +€ jt 3$ 20 20 2e db >< # О аР O 
су iS e. nm OS ч еч Z ++! πο ο DU D A UR DIDA T -- 
e YQ ww eww SSE > «<> > 2 cac О= бахаа х > > 
co Z 1 1 {ΙΙ Ὡς ee A A AO ! — 2 ! ΙΙ C >> 
puc ISI MO C299 COD. +++ кн š H H WH RR HH H ll HH H H IIOU Z + + 
nennen Z m „к=к к= Z Z Zun.0=M оо ми Ом Ὀν ιο ο Z u ZZ 
oO a nr ara en mr js NN NOOOD =< = = = = = nn ern rn rn RO) — OO 
þe << eme ο μι- < -»»»ει og gd» Hoe MM MM ——MCSIOJOJCINIANIONIONEI| UO Ii 
O 2 =— νι OD x > > > = Q. Q. = —4 — — — io od or e pr, por —— C dI а.а. 
O + ПЖ ОСОО) я bBLBELFCOOOOZXIEZEIEXIEXEXXIEXLEIEXExXXOO O OO 
- = ө 
ο <) O O O 
we! e eo f 
WOW ως) WOW 


69 





ου. ουσ ου οι e eee gO C eee oO 
OEM DAO ч СЧС ЧО ЊО С О ety ++ LX O — O Cn O — (NJ G; ШО сос O я см УС CGH OD PO nn 
СО С OO @O G. О Gs" Q." G. Gs O. ОУ О Cs O O O O O @ O O O O me = а СС еее Сее СУ сее се en 
MAY LAY LA LALA LA MA IA LAMA LA IM LA UY OD. XO XO O O ΟΟΟΟΩΟΩΟΩΩΟΟ ΟΩΟΟΟώ ΟΥ ΟϊΟΟΟΟ 9 ῶ 
Ἅῥἑοουουυςθοωοςοοοοοώώσσ ο m οὐ οσο ουσ ο ο ο 
Emu rca κο ο ο ο κ ο EL ο κ ο ο ο κο ο ρα ο ο m aq e. a. ua. «uL aL uu. 
«T «T[ «I «I«T «f «f «T «T «T «Í«T «( «TI «T «T «T «f «T «T «T «T «f «T «T «T «T «Í«f «Í «T «f «T «T «T «T «T «T «T «I «f «T «f «T «f «f «1 «T 
000000 0000000 000000009700 090000000 0009000 010 0000000000 UW 


>> 

>> 4t 3* 

% * > > 

>=> + # 

OO 

COaaa 

Go 

OO e° e 

e ο -- 

DONN 

111.1 

XX >> 

+€ te St +€ 

оссо 
> > > > ООСО 
te 3t 3t 3t O O OoOO O 
>< >< >- >- 9 99ο т 
HAr со со сосо 
C O OO — pad mM HL = ZZ 
SOOO + | + ! Ez == 22. = 
оосо > > > >< Re» mre πι μι ο σι σι Z 
e * * е 3t 3t +£ 4t =Z RNR ZN be U 
AAA OQ >< >< >< >< B me oR ON Z RH OOK NN NR RRR Ne EE 
I l! | ! ἐξ +€ Hr λα. παν ολο πω ον ο ο ο τω <км 
> ><>- >- οσο ο Se ASA OA AO EO AO DO e 
jt 3t 2€ 3* ΟΟΟΟ <> LW KN RRL KOO OW Ee =Z 
ооо! ο οσο UI + eP W>AWF FO PP ZOOIXIHI ZZ 
eS ERN °. е ө е ^, C T7 vt »236 8005 0390. »»ZZ4- 4€ T KR Ze ZI» 
ооо о шшиним ~ С С ¡AAA DIA — — > = ОС >< — 2- Z <=- xç — — r— — 
° ° e eZ2ZZ2ZZÇCGNCGCNI | | C3 X >< *+ — C = > + >< >= O O) = = = Z Z — x Б 
NNNNOOGO! LOS Oe ZKO h tt COOKS ж 
а? ||) чо [| че [к == | Ж єс шшш 915 1 _ 20 mm 4 DU NH MO TI 8% 


DZ EZA rE ооа жо ж Es 
J а оола а.а ОСоа а Oxa oF 
б σσοσοςς-- t шоњо Il (Qu! puru O 
<> € e° o. N l а — x тке gg pe | l HH тп О i 
| d (Sree (amu OH ФО H i —— — H H Ig il ΙΙ UJ 
UJ LJ il wi M ϱ 49 è olj || Hu u 8 @ —— —— — λα. ντα. ete | ο» dm m am ete amm «πρ 
eL XXI» 0 80:00» x » » 0J 0614-09) Qe em СОЎ ОО К O O EEE EN en d N DZ 
ο ο σα σι |! II lH LLJUJUJUJCO (| e "e "e φ. 9. ο. 9. οι φ. M 
На! ZZ ZZ OOoOOoamcoozmiiHdlfee—ee—M Me MM 4 С С СО МСМ СУ СУ СУ СЧ СУ oe 
mm a wd tn dd ed ed ed ed ed ed red went ed ed et wed ed wet Ot ad CO) 
O Ok kk ki ik ik kk = x = = == == == = = = = = == = = = >= = = 


© 
LA 


7 
U 
5 
9 
κ 
x 
x 
x 
() 
0 


H H H H ү 


70 


TRANSFORMATION IN AJ(Es J) 


FORM JACOBIAN OF THI 


C 


NNO Z 
e “(Y = 
c m LLL) n4 
Hn H rs H 
mee = | XS 
CO 
OO =O 
wd ==] = σ--ᾱ 


SIE 
aja О 





QOOOOOOOOOOOOOOOOOOOQOCOO c0 nstuot-0) 0 r O =í n fun O O OOO OOO OG OO 
PFUIVOM ONO ч сыс со ОО = с MAI ГМ O O O O O OD” O a O намес Ом 
EVN COS OVO OO eS S Finn UAUA UNIA UAUN МЛ АЛА АЛУ АЛИДА ИУ UN UA UA UA UAUA UN) U) U.) .O .O -O -O О О хо О 
о.о ооо ою ооо ошооо ошо ошоо О © © © 2 © © о ©‹О О О О ч) О ЯО О О ОО Охо Охо ОООО 
@ O O O O O O O O O @— C, O O O O C, C CO, O 9 9 999909 C OCO O O O O O OCO O CO O OO 9 9959 
л > = > κ ο >Z == z =Z г = > AA aaa a OG 2 ZZ ZZ Z= === 
<( <( <í <( <í <( <í <í <í <( <( <( <í <T <( <( <( <í <( < <T < — = = = = = = = = = = = = = XZ <( < (T < <(<( <f < << << << << 
C)OO QOO QOO OQOQOQUOO[OD OO (QUO ]QOq(QUOoOoOooOooOoooO050500530005500O0UQoOOOOOOOoc( 


FUNCT WRT X AND Y 
BT*D*B 


= 
° ο 
Q = 
ex = 
LU < 
= = 
= e 
--- M ο Q 
— ~ = u. 
Be 2 e LL. E ^) 
» {9 O x =. 
Y ч ë w ~ et 
~ (0% Uu == a 
O O < uy = - 
c O + > m . 
O < ~ red = E 
ο =“ N = ` = 
ο) e < e LU 
+ ш = > bd > 
~ г ~ ч -- a 
x = n 5 a = O 
- < men ш =ч © = 
- uL | Басю O = wy O 
w O ~ Z ONNO < N = ο) 
= N MM. Nam LL > N 
> m 9. — πρ“ ey — C) ο. — — = «-κ., = —. — = 11! 
+ Z N < Naeem = = - =e кюч —— T 
ü f ww eN] a 5.6 e - e o (Y = e ne nn jJ 
— Т ο σ οι н = м м O кс) bk ANNAN IE e CON ου 
- > +} —(<(— Б NZX (Y zo Z XXHXX Z xXXXx o Z 
ww y m~ Z <|  “O% Q< ο. οἱ.) ο P -- С) - = Z — => = ο. == = > ” 
"0 Uj ~ I Ilin ZZ ρω I =— HO m COCO m OOOO v) =з 
Ihr = ~ mm "H H H H HH ze 9 Hx d e lI H H H uw odpm pn nu t Un ή 
T 11) A ^D) ed OAS edt Ol =ч Fd UT) mM — 11! a) ome LL) = | κω num Dre | a — “ο Lit — 
~ O ~ O anran fp 9 0 — JO JD πριν” υπνο Z 
7 230 =“NvV a OO Ὁου DQ *(0OO-zOitiQguuurc OCHHOWWO ш о 
~ UY Ll q www Z Oo or” к 2 ОШ н СМЧ IN im a po € O) O AAA Ll. O 
= FP II — 2222 κ ο» - Z OEE CJ Il ll Be we ew ORO H +... o ч И 
-» q bt pg ΩΙ >x > - O ασ Ον) + O> ANMAN FF 
— _ F e 7575705 Οσο 00 uLIHIOI QOOtJ-—--——-——-—Cogouw.--4—-— uU O 
< > C Qd Aaqa ссе: —4 Q^ O05 qQ O00 c0c00 c0 00 60 Oe — G CO CO О 
ы LW =. = = 
= м) > ac τ e e O о и о 
O < = O DAR) UN (YA UN O O O 
” ο) — U. N u. D N N mou 9 
DON ОСО DOO 000 QUO 


71 





MED DODOOUOIUOOOOOOOOOOOOOOOOOOOCOOOOOOoOoooooooooo 
OO O^ CO mt AI UNO OPEB O O rm CN cn FIN OM 00 ON OD e0J 00 df US «OP 090 O8 O eC0Je00 FIN OM DHOOM 
O SOR Re ММ ММ ВМ co co CO CO cO CO CO CO (Q cO Q O. O° O: Q: O° Q) OOO ON O O O O O O O O O O =— mt em 
DDO O O`AO O O O O O AO. O AO O O O. O O O O.O O O O OO O O O O 09 00 00 oo Poo Poo Piro PP PP PP PSP P 
Μι ε«εοςοοσοσσσσσσσσσσσσσσσσσσσσσσσσσσσσώσσσσσσοσσσσςς 
ЕЕ Р ΓΕΖ Z Z 22 Z 2 ZZ Z Z ZZ ZZ” Z Z LALALA LL LE LLL EL ZA ZA A 
a <( <( <( <( <( <[ <( < >” > Z tt dt tt id dic <( <( <( <( <( <( <í <( <( <( <( <( <( <( <( T L L L L 
QUOOO]Q][UOO(OQOQUOOOOUOOoOO0 ΛΈ ΘΟΕ сакко Q 


+ = + 
+ Бс + 
кс.) < ze + 
K< SK + 
Au. LL) t= + 
Ἡ σος OŒ + 
m U xL 3t 
E < zus + 
+ ΖΞ κ 
+ш чо Ое Nọ 
I WDe Œ e pa 
Hee £ —'uJj Q) <( c OH 
х O OxnxZ-_ Jr m% 
κά οσο) Ιώ 3t 
ioo и) Ou >=% 
ji Οι”. «x 
++ сс. уж «ul > 3 
FIL сро.) + 
+ Oren IO uw at 
3 uuu ur NT" 
Hemi Ler о 4 = Q3t 
+ Œ ke OUJ -- « uJ ОФ 
Фф e EUe mZ = T + 
Κκτουώςμτο KM os 
jx RERXIO => ux 
+ ZU w < шои x 
Fe po et (Dr ας ишн EB Н LLI 3c 
ο iz Z το awn i-i 
+ + Ці MILE <I> <> 
= Home IT P= Rm D O 
ho εαν ο λα -- Ox e 
e че > шюр Oz Zzul «s e 
-- + Fu zu” D WHT SCH e 
> = ки Н исо = = O% м O 
- O ZU =J FO UM b= ἠὲ H 
>< O око udXdxru "uu MX N 
~ = O υγ ο» ο ο. μι WOO O* O 9 
= ο e IX - 7 02.3 Un о. + O = 
Uu» ο. ey ese Um <í <í LLI UJ O. US p= rT st O =з 
< ` M ~H CA ZINJ Ode Is Q м 
nd) --- = tHe Tt OR RWS UDO <> -. 3 
uJ aa -- OH эс Du >+ ο 
+ + == om HU) uj«tz-—uJu» ШШЕ «+ Ow ~ 
-- -- ο. N UENO I 2% su > 
- Y =, ++ ош ос “ο μη eLO + = c 
- ec ~ О =s O = e oz - 27 Om > 
` — ®) ` —- e—dtQ.UUuu»Xu A u. < Zt «x Road < 
~ we D «t SESH >) > HAZ От =% N = e 
— + νέου μοι UJ =+ O C uJc5 
ceo тусу ~ Zaq eto > Zear WOW Ze MM (Sm <T 
> $ — a ч чы МЇ" A — o0 aaa a са e yv i СЕ оса 
-- — ον E DAR OA Ιώ δι ADS UK AZ A ES 
9 = = —) E ee OOR ZZL + OO ON — 
ο» Ec — » CY wa} we 328 UJ — у е) Un 3t e CD жо a Mea 
TOs ZZ O+ = Side 2SYOLN > = www we wet Z NAZ ο 
PY wee e AX ews чә > weil we UJIN e A Z i— C; << >< cx түс [>] _— > 3t wl (em @ € ^ > 
Ueli N «-»-----ς το τας re ea DC nn C) EE ZI c a 0 ~e mat 
κοσμο Hi = αρ Molt mate DIO I шаи < 
ES O SL ee N DE ml AIA (55: DIJO µια + (v <( 
~ = = — < 12 =з = pura ud - δε rr, нб лда. = —4 — Δ΄ «с, с. СЭ w^ 71 
0590979050) «Οσο | «в шр > SOF eK @юш‹—ш®ш ~H It emodtaatones 
О «О *00 eO *Z OOO 6 оу OMWOHD Aw = O DH СОО он κο”, 
ΙΙ." ες» сє +С) а. και 2 < Z= = $ 22 1 2MAYQDAO = 
мо чече d — ANDAR UDI ANT Oe WNMO | os ә 
O-4O0-o000Ozuloozuzozogdsuo «с к ц 2а. 5 c == O— U. LL O U — 
C m CO m CO <( CG) < — £ O) σα σας ος) 2 Uy 0 + UJ ох < < < rJ Z O) O) —— OQ) — < 
3£ U < T £ < Ú u. Z ОФ 
e © © © +. — CO БО) — 623% O 
O O © e 2 jt e 
u^ O © [-- + it e 


(QUOO[O QOO(QOLUOOOOGCOOCQUO 


12 





ννν , οοσσοσοσοοοσοσοσοοοσοοοοσοοοοοοοοοοοοοοοοσσσσσσσσἒο 
GUA DA С О ССО ў О МОМ Ὀσ O OU eo Fin O S O t O SLA OM ORO SAM FIN OM ODO = 
ача а (N (N (NCN CN CN CNCN С\ СМ 00 00) 00 (0 en e e eO eo eo X. JF ST SR GS ULA LA LA LALALA LA LUN O O 
М М P P Pee Se P P Pe Pe P P Pe e Pe P Pe P М М See P Pe Pe Pe Pe Pe P P М Pe ee fS IS P SS IL IS 
NEDOOOOOOOoOooOoooooooooooooooooooooooooooooooooooooc 
A 7 22 2 2 Z A DEL >=“ ώς ZELL LEI ZELL DZ PP ZZTZ 
<( <( <( <í <( <í <( <í <( <í <í <T <í <T <( <[ <í <T <( <í <í <[ <í <í <T <í <í <í <( <í <í <T <í <í <T <í <( <f <í <í <í <£ <í <í <T <T < <1 
UU OO UUU UU UUUUUU UU U0 UU OUUDUU UU UUU UUUUUUUUUOU UU 


GO TQ 300 


I 
2 


= oO 
м OY » 
ma mH NI 
ul Z— e! 
m ES Y > It LW 

pe οἱ!) ı —D 
OD ew IC + Z — rL << N 
O TKI Ne 6} DRO STN 


pues |] ме 


AUL s? NL) -AAXACIS,.K) 


NN 


"— Кеш a 


‚ 51000) 
5 
5 
5 


έ 


II —J3 Z—=—+4 X LL] = O mr Q. >` +£ 


FOLLOWED BY A 
AFTER THE 


(B), 


DECOMPOSED BY THE SUBROUTINE LDLDT 


Ἂς ἃς ας ας ας е с д д ας с о д ж ж atate ote ote ate ole ole ate otc ote oc otc ЖЖЖ Де ЖЖЖЖ ж 


ORMS A FORWARD SUBSTITUTION; 


x 
it 
<I 
> 
uy 
z — 
(ἢ 27 
eo > 
«τ Uu» 
= LLI 
`+ E 
m > 
— + 
Nim α 
lo C 
Ce τα 
h Q. 
= T++*+ OZ 
= | =+ LO 


= м 


DODUILZLOJ-OUT-O0OD Ζω, ὠνα 
OQ. << C сәс оци ” WNN Y mO) e 


e NY 
SORO ОСО 
EU C» OO 
Ce) f OO 
in 


“г< < 
3° — cO = 


, 


«1. J. 
a 


ONE CALL OF SLV PER LOAD VECTOR IS REQUIRED 


CALL TO LDBW 


m 
1 


AFTER ON 


EEE EL EEE 


972 
ж жж 


zt о 
С 
e CD 
«x + 
үк се pang 
ho =>) «- 
Y ο. - 
С) ч E 
= (D W— ~ 
e «t ec] a A 
E Ole] N + 
bond Th ! II Z ~ I = l сс 


CD th — “Зз є ЈОД) ева Јн et 


ш 
Ir 


l NM 


a — μα N < de — 
>% Око | 002 ОСО 1 
СОФ T OSON eI -OlI o 1 И 

Je Zoe 00M M LT OI--O Za 
Cie ii N ee є— μος 
Wk ZOHO I LL O O— Oc — 
Q Z OOO 2-030030 0-600 


Qit 

C OTO 
ic O 50 
E | су 


Суусу OUO OC UOOUOUUUUDOU 


13 





ΟΟΟΟΟΟΟΟΟΟΟΟΟΙΩΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΩΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟ 


Q eC r SL .O IS O Q* O —CN CN rO Ωω ο COS QD чтим 


eA tm DANO FIN OM ODO HOIM FINO 


шаі ol ul eo e cd cet mm OSOS OSO С\! С\! СМ СМ С\/ СМ СМ СС СС) СС) СС СС СС СС СС ССС О Р ОО ОО EE MMMM IN 
сз со со со со со со CO OO OO CO CO CO CO CO CO GO CO OO CO CO OO СО со со со со со CO со CO CO CO CO CO со СО с ОО ОЭ со СО с O OO со СО со 
@ O O O O O O O CO C O O O O O C, O OO OCO OO O Ce Ce CO CO O CO O O OCS O O OO CO O O CC O O GO CO OCO O 
уг уг г суг Сг eee ze ALLE EC СА СЕ сЕ СЕСЕ СЕ СРЕ 5 5 ШЕСЕ EC EA СР 
a< <í <í <í <í <í <í <í Z <1<1<1<1 52 <1<1<1<1<1<1<1<1<1<1<1<1<1< < <1<1<1«1<1<1<1<1<1 <<< <1<1<1<1<1<1<1 
EE ee ο ο SSS = τω 


осо 


cm 
ooo 
NOO 
NOO 
m 


© 
O 
e 
a 
NO 
5 
= OO — 
cN CJ Cy GO -- 
7979 оо nO 
о wn е ο O Juno 
<< NN e © 
(UN a мо uñ 
-»-- e mm t w «f 
w EE wem = 
ELE zz OC Z— © 
ZZ e OO ~ «(4 ο kt 
эч nun MAN δ mo 
oz NN mei) Yew CO 
+ + i= c ku ~ О 
— MÀ δι nn 2255 о ~ 
ei AO ne Ne ο O 
со MA NNMNZZMO O 
{= Гү μι ο) κιν ϱ =ч 
— — = f ww + | >X<r(<— O O 
zd > Be) =» ° 
px --Ὕμρ. πη Ee n NE 
d NZ m~ 077 πλου ° 
LU NAZ mm wZz=nmanann e Шш 
O. NUDOS * и са = x yx; — > dc 
ZN Sr il DO ZINN I mu xD N eZ 
rie (JJ. AA O Al N SET 
[TAT em en TI 


11 νι | 
ο) ο — < < << 


- ο (Y |] par 1] —> > — > = j+ Lx vO << 


s ry Z Ç e rn ¿un P 


NI 


© 

e 

e 

© 

N 

м 

оо — 

OO — 

оо un 

σας 88 Y 

rt (UN N 

a — = 

3t en чь? 

O mma) [= 

O cc. << = 

© ~ «(5 э 

e MNK Z 

N =e Y 

„ Бо a 
OO > 525 =P 
OO 2202200) ο 
e e ИО аан 
NO + I >x << —— 
= + > we IE 


< > = ry (rC Q O. Z 5 
CD «X (Om >< >< (/) Ç) 


ας 
m 
ον 
= 
сан? 
= 
< 
э 
Un 
ç) 
e 
e, "πρ 
ne Ὃ 
Cs ο 
e wd 
I! — T 
ngt = 
ο. е бм 
O E 
— 4 m) 
eo ο 
po po »-- 
ч” 0» u 
>< > 
ZA = 
~ = 79 
v. = 
N е Y) 
€ 1l y 
=< > 
(N =ч 
= νο) H 


Y Pr ee RY Ja Oe OO 
o jA  Ι Ἡὶ oZonozZo 


--- жоу 
С^. μα ο nn 


Dd erde my 


πλω το “τ, πο 


— 


ο ο ο ο Κω ο 2224 mm em ο |) СУ |! 


mM Z= ә меше о осет ЛЛ СЕ сту сус, ` Jg ΜΓ О) 7 ССОУ ССу СС) СС et N ον φ 
О 449407) Z Ы әла) и) В a ee Ne) 2—20 d 
OVO www IONOV i I O~ CDC < Ом н Il Omer meo e e O) se jua o O) e 


OOF I III zolomcuos» Il meet) LOOOMU> > i kiki UlCOUI «Ul OL 
<и себ исә 22220 +22 ноо 220 G; II = II ZI DDr ZZ ZZ Zr Ol. it. 
Z —— 0.2 Z nne ZZ Z 0.2 Z "5-5". οὖ» »- 


Re < = 


с) Осо ZÉ Z Uu. O (O O (= <( >< >x — (n (y (O AL, TO <Ç < x Am Z Z Z Z K O - G K O < 
OLLU ee C) CJ Y) UN UN CU C930 CoU U^ <t () O — ч í (93 (5 (O (O < A n n m = G; xu. = G = 


e 


CC 
(A72 
NO 


© 
© 
+ 


το. 


OO O 
nun e 
++ ιω 


Cc O 
оо © 
OO = 
ex 





@ O O O O = (NI (9 + try O I — O О чм омот O (ju rn O OO OO O O O O O O O OC CO CO CO OOO OO OO O 
P< GO O> O erm — — — — —4— —4 —4 (NI (NI CN. (NI (N QUI CUR CO OO CO SELLA DH DO SAM Ωω OM ODO SAM EIN 
UNUNUA JSS OO «OO OO ODDO DO-QOO OD OO DD Od O0 ШОМ М F F = F ММ М cO cO OO C co CO 
QO cO cO cO cO cO CO cO cO CO CO € CO cO CO cO cO CO cO cO CO cO CO CO cO OO CO CO CO CO CO CO CO 02 mO cO CO cO cO CO cO cO CO cO cO cO cO cO 
€ O O O C O O @ O O C CG O C O CO C CO O Co CO OOO O C; C O O CO OO © O CO CO OCO O OCO OO OOO OC OO CO OO O oO 
zz z Z ZC CO OCyYC C C AC CO C C C y CA C C CO CI CC CY C C Z Z Z Z Z Z ZZ Z zz nz ern 
<í <í <( <r <r > = 2” > > >° >: = Z= >x > > > >: > >` = = X> > X <( <( <T << <( <( <í <T <T <( <í <( <T( <( <í <T <( <( <T( «f «fp f 
QOO OOO000000000703050030Q00€30 O000 0O0O0 OC OOOOOUCDOOD O0 DO DO O VDD 
>< 


>< = Un» 
ect < e ш 
> > = - c 
ne me x e > 
Њи) “= LL en > v 
== Ne ο. + = 
„ч a j a e = Ф 
ο нм >» = O N 

“ Oes ш Z od e 

e a T ame = ~ 
ш «T > © W 
το = ο _— N 
ке WSs ο > о ч 

= τα — м (C 

e = YN x~ = о 
h 9 = т NN e = 
ex & e er ~ ο. 
m EX >< N o + = 

e δω el^ e e x+ 
Ne m N - ο ο“ O w 
| ` w “м, Ne e nj E e 
wu» e = Y) -— N 
| œ | >e = x = <t z o 
nz = i ZZ O m 
f ejN) е O P Q ~ 
NINI Q. C mo Lon] O< - Ua μι 
le u tuu N “905 w Oo C 
uj eO «fe C μια - “eZ 0 O 
| x<— | u I e Ui e — (5 = Y 
eC orm Axa Tx Tx AID = e 
| eu | οι 2 Tu ecu = 
he > «+» = = = e “- ~ «OZ © 
t | œ < = ~ Le ++ mi or OL ad m 
ux uU e << e QN + Z= τ- 

65 >< = ¿"J od > od wen) 6.27 e 

Oe ο. = LLJ II nu H C FO = ° 
LL] = < ЦО) >< == M5 E т 5 «1 -— «I» Ὃ 
| GN | G. N OO Ni e j= e Ос ос e 
Me — DU) 4 NN Lulu) — < =ч „= e 
| == |= = N We S qu) ә 200402 O 
«Тусе < es t atm = = 9. eQ. e ej Zz о LA 
| — | x + Ux bd ‚ш = еъ <] eC) — 
OC e» Oc ex ο ο + === = x м СОЈ tb 1 e N 
le «| e ” e ule = ue Eb Ol e TE © m 
We ule c Sur Oe = Ox zZ zo +Zz + і 
IZ e |> e band band <{ “3. <P 7 {οι = e 
>ч > 2200 *== wa ce 2, ny e Z DIN = =— e O 
(Os jm» = шз Ne Ue y σι SENO o = 
1-5 οσο. az > O N >Q -- => TOMO m 

< эх EB -- <í > --- «Ὁ e Qw Wd zs -- © 

ee Po m ch Noa - - med oe Dem гоо eeN C) e O00 
e = dem > minu nn Ty .- en =~ (с Z PB = rr O Te) 
>x > s> = = es =l || — >< >< Lem m I U Zn €D κο ο UN E oom (N co ~ 
„ЧІ wm —4>< = sitet — CO = UN ο) Ωω) ο) nous σσ οο Zu N HM cfe 
= aZ anz loo lo ФО = = FOO e —ə ШІ) ММК Сә erii 9. ο jl ο 
ММ еа фә (VM MO eO NX lO ULA JEND ο) O AS 
“ШЭ 5 “54 І >" 6 6 МММ осом он Z (0Ο) 2ου "πο νιον N 
an e. = vermis! эм N Гу м e єз ч ОЎ наве ч М. а мн С) а н оО) ч і 
w V~ ш OoZom—I209- MC M) O = C.O pr AAA PAR OR Oe OI KIN e 
> =» 9-LZOl——-—-——— u Cohen Dun Жш = с-ф к NE UL eu — 
«t1 ec A A Ou. Duo! e ZOZ+H Ome OOO uD 
> XX Очен Бо сое И SUCA 2 Modo So = Cerny H Ων li ο δρ» 
x ми Ne Z 2X — x - o. мса = 50 но БООН X= к= 
Oo eC HULOWONZOCKO »5oO0agoaouoncauzozooozs-ulo""uuoaoguicoiuooocd 
U. e {} e CQNOQONNOZIL «Сн Ж. СУ 0CUJUO QC * LO) tin x 09 mm CO — C CD — C5 

AN N eU і ~ 
O O O оо о OC © O Onl оо 
O O O AN © OO O O ind cN cO 
e o u^ nn y CC) O O — rn 
-- en JT + T О М D Ον 


76 





τ ος ου ου SOS SOS O OSOS OOOO 
OR СФО сыс, ИО МОО О О AAO FID ο 00 PO AN ASA οσο O E LN O Po OO O Nm) 
DOIDOINNNTNNNNNNOOOOOO OO 9 Am ee rd OI CN СМ СМ СМ СМ СМ СО СО) СС) су 
CO CO cO CO CO cO cO CO cO CO CO COCO со О\ ОУ О\ ОМ ОХО ον σν Ον Ον Ον σν Ον σν σ» Ον σν σν Ον Ὃν ө 0 0 0 е Тө Хө 0 Те о Тео Тете 
OOOOOOOOOOOOOoooOoOoooooooooooooooooooooooooooooooo 
ES uuu uae ME E A Z= E LAA AL A AAA LA LALALA 
«L«I«pI«(«p «(«r «I «I «C «f «f «I «f «f «I «r «f «I «f «I «f «[ «T «I «c «f «f «I «f «f «I «f «T «f «T«f «I «p «I «f «I «f «f «f «f «c «1 
SOODOOOOODOO OOOOOO OO OD OO ODDOO OO OOOOOOOOOO OD OO OO OOOO 


-3 
«t 
VY od 
< 
ео 
«tx 
N 
х < 
om /) 
mitt 
oy 
б. ~ 
С? =ч 
uJ e 
a ts EY 
~ LLJ 
y at 
Dr 
DNS 
+ © 
co 
{ 
Qe 
9 v 
QN e 
AQ. 
С 
> «f UJ 
x > Mm. 
O Xx tM 
> # O - x 
>< -- % z a AI 
O — m~ «T«r inl end va 
+ e 7 NN «t «t K 
== PUE e Jt 3t Uu» on + 
= O O — 34 46 A 
e cc ш => Qe «=ч e) 
ο -- 4 O ο orn CP <I 
ш ----- e = = m A шо 
»- Y Y © OZ (Y mt + 
~ O 0 e сс СО Oa. w — 
Y MN ο) „~ СС h CY wa `Z <t 
со | +I со «--ν- шх On 
a +i Ον = Ne Du YA m t! 
! SI D ч м со хк Ome CD + J 
=ч = OOO =ч pr 00 эсс << 
= Ao ; а. + со ~D | IO 
ш Om = С) + Ο.. | CIM WJ O? 
— eu I к =ч pr L al LLJ + 4 + 1 
w E (η y IR N alt 
~ о С) очы чо — ~ e O O*o * «r3 IO = 
Ст = = Ρο O = == pæ m Came te at QR Oe Or Owe O 
e OOD - 165». qe ο —y << < EE TX H de c 
= JAN = и m |- > "TT тю m= TITS H# *— 2 Hl⁄ Il = ~ 4 СО < 
—» << == ZI n 5 = wA ZN ιο ον nm Ot NT tT ț 
к=з Il + ιν. je 9 — — ο. | ο. d» а г en ка e ell flo 
w οσο δι οι ο u w (ITAM) "--ΌΟ 5 I QO DO mm 
QOo-—Qg-—-11 1-000 WOam Aan COCO + | JA ο ο Ην € guide ο 
= < + НЭ СУС ex ϱςο sC στ «гч UU Mery 9 6-4 н > πι C. Lla € O LLS 
HOOD ulurz-.)j-—2-«ftNO- Ca —C 2 ш =D.» ED mm Cry CZ > 


О aJU uy OA O) rr OD ir DA ым чм O I| CO O — O = — C *Z coo ule x ululzwxxsxculululzz 
i COCO L и I! ZC COCCO COO I OO li OS OOO и со н FL CO CD mee СО со а COCD i = — — c 
Ib lili Loser Je ee l! i= <í Il «III H I O — 9 СО СО ч — мес мБ сон нн a S 
е“ >С) DO-I-0 Zw > ελ œ AJJI GC II H I `Z σοι η ων σι * x * Z= 
OQ XU JOO YOL JU. JU. OX 0xX«OQO OO. .J« «t uJOdo Xu COLL σσ «20 0) ιο (2ϱ οσο! 
(50) AL O C) TL YH) A DO C3CD CD «T6900 9 Ce 70 et «Tf cO CO «X «T 00 cO (0 «1 «x co cO cO (C5 XL 


N ONS + ο о O N oO ra e O 
M rtr + y n ο о ~ r= со O 
-- ~ = ==] = = == == = = = ο 


del 





mooooodoouoooooooOoooOoodocOdoococoqOoocoOcc0Oocooooooooooooooo 
SUS OPM OD OO MAL EOE OM 00 0 O eA CU001090. 00 0e CD eu NIT LN SOM CO PRO ών οσο ο 
CACO 00 eo c0 eo + + + + Sp SP SR SR S$ ο “οι ωισιωιωισιωιωω ο ο D O .Q -O .Q -.O O OUT =I&Ñ = = === оо 
ΟΟΟΟ ΟΟΥΟΟΟΟΟΟΟΟΟΙΟΟΟΟΟΥΟΟΟΥΟΟΥΟΟΟΟΟΟΥΟΟΟΟΟΟΥΟΟΟΟΟΟΟΟ ΟΟΟΟΟ 
wr OOOO OOS OD DOO 0O0900 OOD OOOO OO COOOODOOOOC OOOO OO00OoO 
AA Z LA ZL ZOZ 27 2 Z Z Z Z Z Z Z Z Z Z Z Z Z σα στα σι ισα σσ σώσω” 
< <( <( <( <í <( <í Z” <í <í <( <í <í <( <T <í <( <í <í <í <( <í <( <( <( <( <( <í <í <í <í <í <í <í <( <í <í <í <í <I < < <I <L <1 <T < «1 
lU u G Gss a sss Se 


CON(1825,14) ;NBC(217,2) ;NSTRE 


:ELCON(10,4) ; TITLE(10) 


ш 
e 
Q 
DZ x 
Am € - 
ecl) E 
> DZ O 
ed om CL »-- y Pro 
AIP Zk = = 
wi ow 4 es 
OSO = = e 
«I» c O O 
оао 2 ~ mij e 
UTC A e = O 
αμ δω ο O ο & 
eC) == < un ao e 
ma ~ aQ a U mo «t со 
NAN "re JH c e 3$ (UN = 
| =—<—s a O ~ MD ~ 
Om + = Z = €. + GOG --- i 
am M Z e a =ч LL) а. 
TO * .-+-0 » ο. — OO + 
| — + — (N >=—— O = th «=ч 
cu Cte ur Το’, μπε Ar AL --- — a - ο 
amar Pe ee О з “ә OO ο -- e N» = 
ο ο tmb I CUNO e - (505 e e ο. a wean ct 
N F OX MUSNA со (=) Su = = к= |-- “<I <I.) 
чыш) ш — > < INO == cN ιν = “η D €) к. ШТ 
OINOZ eu ЧО. с) = na ον к: =a m — wh Oe 
GANAS ZN mm eo - Су а — ~ e — OIJ 
Wer 0. „Jene (OQA οί OO c ΘΟ C C O O qaq UJ 
Z IO VMTN Ne OD tt Oe WOON CN οσο στο το ALOQ ~w wmm 
--ε]--ἰἰ. YY A o ο wt et ος) C — + =y + +O - - ο )-ι OCH Om Deda 7 — 
FINN SNe Or DO ew SG άσσο e wet md me ed a Јес а, 
Dvzzz JznzHoODZ1O ш mm -- ο ο OQ N «(o Οἱ, ноло 
ОнОоо =C) Z > sC Z — II < lI =+ LOO Ol ONONINO οι "oINgHCsooOos- 
AAJESN OZ IA MO | ες οί 1 AOL ee ee ied A O) A O) E LA Re < (CO T Il ο 
COS = = = >= DI О Ос) OI AAN ~e~ lI "= || - СО HOO me 0 —D a Jl 
Z >= CG)O O Z O —Cz+—( UDO κο O LL K L.—:, LL. LL O £ OO O z O O Q Q. Q £ 9- O O X Q. O OQ. O, < < — 
WY mI) OWA к OOS ILD OM = 7200-00 0(00 MO Di > O CD — OQ) (5 << CO C) <( 
2С : 
O N © С» ορ ο (a) O N © 
— =з N M (ҮК ιο O O M 
A qd ri nl = — wi pem e=] ---4 


78 





ооо ооо о ооо ооо ее 
NOIA O A DADO O -—cC0) M uo 
@ OO CO COCO CO OO OQ CO ON GX OY G GY O: ON Os 
ΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΟΥΟ 
eleleleleieielelelsieleieleielel®) 
Μι: = Z Z Z OG Z CO Z= 
<( <( <( <( <( <( <( <( <( Z ачса 
WOU VOLO UUUOOO00000 


m, unm, 


"Ny 

9» » 

um pend 

u s 

T 

OW 

<< 

шш 

οσα 

++ 

rn ч 

N 

22 >» 
Qo ao 
coco ο 
<< zx 

+ + t ee 
+ OO Fr.» 

< < tC 

.J) WO AS OO «<< 
E os < 1 его) 
Uu) e | 1 COO mo oc 
+ QQA ο tu l j O. H H 
C$ шуш << j σου A un 
LLJ =d NINE E a ө Әј ә ὃν 
— + — Ja aod 110.9, 
Ra το σσ tud M 072) 73v 77) 


0 WNNWW=HDREIS == 

i CO o o Sh be О С 
Qa | OOOO | OUMOQUD 
ПСИ поо << сад <4 «ito 
Ima AIDOLTA AX UI 


O O 
O O 
N m 


19 





LIST OF REFERENCES 


Porter, F. P., "The Range and Severity of Torsional 
Vibration in Diesel Engines," Transactions of 


The American Society of Mechanical Engineers, 
vol. 50, APM-50-8 , p. 46, 1928. 


Cantin, J. G., PLISOP, an internally developed FORTRAN 
deck, Naval Postgraduate School, Monterey, March 
w972. 


Software Licencees Limited, TORT 2, a FORTRAN program 
for the stress analysis of circular sectioned 
solids with varying diameter, Swansea, U.K. 


Zienkiewicz, 0. C., The Finite Element Method in 


Hne neering o erence disc a S 


Adamek, J. R., An Automatic Mesh Generator Using Two 
and Three-Dimensional Isoparametric Finite Ele- 


ments, Master's Thesis, Naval Postgraduate School, 
ШОП егеу 1973; 





Allison, I. M., "The Elastic Stress Concentration 
Factors in Shouldered Shafts, Part III: Shafts 
Subjected to Axial Load," Aeronautical Quarterly, 
ον ΙΙ ΠΉΠ,, μεν 1352. 


80 





INITIAL DISTRIBUTION LIST 


- 


Defense Documentation Center 
Cameron Station 
Alexandria, Virginia 22314 


Library, Code 0212 
Naval Postgraduate School 
Monterey, California 93940 


Department Chairman, Code 59Nn 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93940 


Professor R. E. Newton, Code 59Ne 
Department of Mechanical Engineering 
Naval Postgraduate School 

Monterey, California 93940 


LCDR Dennis M. Doyle, USN 


213 Canterbury Drive 
Wallingford, Pennsylvania 19086 


8l 


No. Copies 



























Г TIT ЇЇ | d 


“DUDLEY KNOX UBRARY 











Uu 








